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[i]  This  report  describes  thermospheric  composition  and  particle  precipitation 
changes  that  occurred  during  the  period  of  the  great  geomagnetic  storm  of  20-21 
November  2003,  an  event  that  was  associated  with  the  passage  of  a  magnetic  cloud  past 
the  Earth.  These  changes  are  compared  to  those  observed  during  geomagnetic  activity 
on  17  November  2003  and  during  the  intervening  quieter  period.  The  data  used  are 
obtained  from  (1)  ground-based  magnetometers,  an  imaging  riometer,  a  scanning 
Doppler  imaging  Fabry-Perot,  and  photometers  from  stations  in  Alaska,  (2)  photometers 
from  Canadian  sites,  (3)  NOAA  POES  and  DMSP  particle  sensors,  and  (4)  the 
TIMED  Global  Ultraviolet  Imager  far  UV  sensor.  The  composition  changes  associated 
with  the  input  of  auroral  particle  and  Joule  energy  showed  larger  depletions  in  atomic 
oxygen  on  20  November  than  on  the  other  nights  and  greater  changes  than  are  seen 
in  the  Naval  Research  Laboratory  Mass  Spectrometer  and  Incoherent  Scatter 
(NRLMSIS)  model  atmosphere.  NRLMSIS  does  better  in  reproducing  the  changes 
during  the  great  magnetic  storm  with  its  long  duration  auroral  energy  input  than 
during  the  shorter  time  duration  geomagnetic  activity  that  occurred  on  17  November. 

During  the  nights  with  the  largest  changes  in  composition  the  input  of  Joule  energy 
dominates  over  auroral  particle  energy.  It  is  shown  that  the  particle  energy 
distributions  associated  with  the  20-21  November  storm  in  the  period  around  and 
after  the  passage  of  the  magnetic  cloud  had  lower  average  energies  and  were 
enhanced  at  energies  below  0.1  keV  than  those  that  caused  auroral  displays  on  the 
preceding  days. 

Citation:  Hecht,  J.  H„  et  al.  (2008),  Satellite  and  ground-based  observations  of  auroral  energy  deposition  and  the  effects  on 
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1.  Introduction 

[2]  Geomagnetic  auroral  superstorms  are  rare  events  with 
only  a  few  occurring  during  each  11  year  solar  cycle 
[Val lance  Jones ,  1992].  Because  of  their  rarity  and  the  large 
amount  of  solar  energy  that  moves  to  the  Earth  from  the  Sun 
there  is  great  interest  in  understanding  the  nature  of  these 
storms  from  birth  to  death  throughout  the  entire  Sun-Earth 
connection  region.  With  each  successive  solar  cycle  more 
data  on  such  storms  are  acquired  as  an  increasing  array  of 
ground-  and  space-based  assets  are  deployed.  During  this 
solar  cycle,  number  23,  several  large  geomagnetic  storms 
have  already  occurred.  Two  of  these  are  within  1  month  of 
each  other  in  October  and  November  2003,  and  they  have 
already  been  the  subject  of  numerous  studies,  taking  advan¬ 
tage  of  these  observational  assets,  which  were  reported  in  a 
special  issue  of  JGR  [Gopalswamy  et  al. ,  2005a]. 

[3]  In  Situ  Mass  Spectrometer  satellite  measurements 
have  historically  been  the  largest  database  for  information 
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about  composition  changes  in  the  thermosphere  during 
geomagnetic  activity.  In  particular,  the  Atmospheric 
Explorer  series  sampled  down  to  135  km  altitude  and 
data  from  this  series  of  measurements  form  an  important 
database,  which  is  incorporated  into  the  Mass  Spectrometer 
and  Incoherent  Scatter  (MSIS)  series  of  model  atmospheres 
[ Hedin ,  1983].  These  and  similar  satellite  data  have 
been  used  also  for  many  individual  studies  during  both 
geomagnetically  quiet  and  disturbed  periods  [e.g.,  see  Nier 
et  al ,  1976;  Kayser  and  Potter ,  1976;  Prolss ,  1980].  These 
data,  however,  only  provide  a  snapshot  of  the  effects  of 
auroral  energy  input  at  any  one  location  and  these  data  have 
not  been  available  for  many  years.  Recent  studies  have 
involved  satellite  remote  sensing  of  far  ultraviolet  (FUV) 
emissions,  which  to  date  have  mostly  been  confined  to 
dayglow  observations  [e.g.,  see  Nicholas  et  al. ,  1997; 
Strickland  et  al .,  1999,  2001,  2004;  Meier  et  al. ,  2005; 
Crowley  et  al. ,  2006]. 

[4]  However,  over  the  past  20  years  there  have  been  a 
number  of  studies  which  have  reported  the  effects  at  night 
of  auroral  energy  deposition  at  well-instrumented  ground- 
based  sites  such  as  at  Poker  Flat,  Alaska  and  Kangerlus- 
suaq,  Greenland  [ Christensen  et  al ,  1997;  Gattinger  et  al , 
1991;  Hecht  et  al ,  1989,  1991,  1995,  1999,  2000,  2006; 
Niciejewski  et  al,  1989;  Vallance  Jones  et  al,  1987].  These 
studies  have  shown  that  often  there  are  large  changes  in 
atmospheric  composition  with  only  modest  inputs  of  auroral 
energy.  There  has  been  only  one  such  study  of  a  geomag¬ 
netic  superstorm,  which  occurred  on  8-10  February  1986 
[Hecht  et  al,  1991].  During  that  storm  there  were  much 
larger  changes  in  the  ratio  of  the  density  of  atomic  oxygen 
[O]  to  the  density  of  molecular  nitrogen  [N2]  than  were 
predicted  by  the  MSIS  models.  At  that  time,  however,  there 
was  only  one  instrumented  ground-based  site  (Kangerlus- 
suaq)  whose  data  were  analyzed  and  no  attempt  was  made 
to  relate  those  changes  to  auroral  heating  rates  (local  or 
global)  or  to  global  auroral  morphology.  A  recent  study 
suggests  there  is  a  relationship  between  local  auroral  heating 
rates  and  the  [0/N2]  ratio,  but  the  data  are  still  too  sparse 
(especially  during  large  storms)  to  determine  if  the  relation¬ 
ship  is  linear  or  of  higher  order  [Hecht  et  al,  2006]. 

[5]  The  20  November  2003  geomagnetic  superstorm  is 
particularly  interesting  as  it  was  reported  to  be  not  only  the 
largest  geomagnetic  storm  to  date  in  cycle  23  with  a  Dst 
index  of  -472  nT,  but  it  was  also  associated  with 
the  passage  of  a  large  magnetic  cloud  past  the  Earth 
[Gopalswamy  et  al,  2005b].  Originally  defined  empirically 
in  terms  of  in  situ  spacecraft  measurements,  magnetic 
clouds  are  observed  as  large  transient  structures  in  the  solar 
wind  having  smoothly  rotating  and  enhanced  magnetic  field 
vectors  with  low  magnetic  variances,  decelerating  speed 
profiles,  and  low  plasma  temperatures  over  intervals  of 
many  hours  to  several  days  duration  [e.g.,  Burlaga,  1991]. 
Such  periods  are  easily  distinguished  from  the  ambient  solar 
wind.  The  simplest  interpretation  of  the  magnetic  cloud 
structure  is  a  cylindrical  magnetic  flux  rope  [Mulligan  and 
Russell,  2001]  that  is  either  expelled  from  the  Sun  or  formed 
during  the  process  of  a  coronal  mass  ejection  (CME).  In  situ 
observations  of  magnetic  clouds  can  often  be  related  to 
distinct  solar  ejecta  [e.g.,  Gopalswamy  et  al,  1998;  Webb  et 
al,  2000].  In  particular,  magnetic  clouds  observed  at  the 


Earth  are  most  notably  associated  with  Earth-directed  or 
“halo”  CMEs  seen  in  white  light  coronagraph  images 
[Gosling,  1997;  Zhao  and  Webb,  2003]. 

[6]  There  were  many  ground-based  photographic  observa¬ 
tions  showing  significant  red  emission  (e.g.,  see  http://science. 
nasa.gov/spaceweather/aurora/gallery_0 1  nov03_page2. 
html).  Such  observations  suggest  the  presence  of  an  unusual 
distribution  of  low-energy  electrons  well  below  1  keV  that  are 
the  suspected  cause  of  the  rare  Type  A  red  auroral  displays 
[Robinson  etal.,  \9%5,  Sivjee  and  Shen,  1997].  An  earlier  study 
of  optical  emissions  during  the  1 995  and  1 997  magnetic  cloud 
events  showed  they  were  in  fact  associated  with  the  presence  of 
low  average  energy  electrons  that  are  one  of  the  hallmarks  of 
type  A  red  aurora  [Sivjee  and  Shen,  1 997;  Steele  et  al ,  1 998]. 
However,  those  studies  were  confined  to  a  single  site  and  only 
inferred  average  energies  with  no  information  on  the  energy 
distribution. 

[7]  This  paper  addresses  several  questions  associated 
with  composition  change  and  auroral  precipitation  during 
the  period  around  the  20  November  2003  geomagnetic 
superstorm.  (1)  Is  there  evidence  for  lower-energy  electrons 
or  other  differences  in  the  energy  distribution  during  the 
20  November  superstorm  as  compared  to  the  precipitation 
during  geomagnetic  substorms  that  occurred  before  and 
after  20  November?  (2)  Is  there  a  connection  between  the 
passage  of  the  magnetic  cloud  during  this  storm  and  the 
low-energy  particle  precipitation  as  suggested  by  Sivjee  and 
Shen  [1997]  from  earlier  observations?  (3)  To  what  extent 
do  the  0/N2  model  (Naval  Research  Laboratory  Mass 
Spectrometer  and  Incoherent  Scatter  (NRLMSIS))  predic¬ 
tions  and  actual  observations  agree  during  the  period  before 
and  during  the  superstorm?  (4)  In  particular,  do  we  see  any 
difference  in  0/N2  between  the  model  and  observations 
during  the  time  period  of  the  magnetic  cloud  event?  (5)  Do 
any  of  the  observed  atmospheric  composition  changes 
correspond  to  inputs  of  auroral  Joule  or  particle  energy  into 
the  thermosphere? 

[s]  To  address  these  questions,  this  study  will  use 
ground-based  photometric,  magnetometer,  Fabry-Perot, 
and  riometer  observations  from  the  Alaska  chain  in  Poker 
Flat  and  Fort  Yukon,  photometric  observations  from  the 
Canadian  chain  at  Pinnawa  and  Gillam,  Canada,  particle 
observations  from  the  Defense  Meteorological  Satellite 
Program  (DMSP)  J/4  and  J/5  sensors  aboard  FI 5  and 
FI 6,  particle  observations  from  the  Polar  Orbiting  Environ¬ 
mental  Satellite  (POES)  Total  Energy  Detector  (TED) 
[Evans  and  Greer,  2004]  aboard  the  National  Ocean  and 
Atmosphere  Administration  (NOAA)  16  and  17  satellites, 
and  FUV  observations  made  by  Global  Ultraviolet 
Imager  (GUVI)  [Christensen  et  al,  2003]  aboard  the 
Thermosphere,  Ionosphere,  Mesosphere,  Energetics  and 
Dynamics  (TIMED)  satellite.  These  data  will  be  used  to 
study  the  period  from  17  to  21  November  2003,  which  is 
just  before  and  during  a  geomagnetic  superstorm. 


2.  Experimental  Technique 

[9]  The  approach  taken  in  this  study  will  be  to  determine 
composition  change  primarily  from  ground-based  photo¬ 
metric  observations  from  Alaska  analyzed  with  techniques 
that  have  been  used  for  many  years  [Hecht  et  al,  2006].  The 
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Table  1.  Geomagnetic3  and  Geographic  Coordinates  of  Ground-Based  Locations 


Site 

Geographic  Lat 

Geographic  Long 

Geomagnetic  Lat 

Geomagnetic  Long 

Fort  Yukon 

66.6 

214.7 

67.3 

-94.7 

Poker  Flat 

65.1 

212.5 

65.44 

-95.66 

Gillam 

56.37 

265.40 

64.54 

335.14 

Pinnawa 

50.15 

264.12 

58.48 

334.19 

aThe  magnetic  coordinates  for  Poker  and  Fort  Yukon  are  in  the  APEX  system  [Richmond,  1995].  The  magnetic  coordinates 
for  Gillam  and  Pinnawa  are  in  the  eccentric  dipole  field  line  (EDFL)  system. 


average  energy  of  the  precipitating  particles  will  be  obtained 
from  Alaska  and  Canadian  ground-based  photometric  data. 
To  obtain  both  quantities,  typically  two  auroral  emission 
ratios  are  measured,  each  sensitive  to  composition  and 
average  energy.  Hecht  et  al.  [2006]  suggested  a  variety  of 
auroral  emission  measurements  could  be  used  to  obtain 
these  quantities.  For  the  analysis  in  this  study  a  number  of 
these  techniques  were  used.  However,  during  the  analysis  it 
became  clear  that  most  techniques  did  not  provide  a  good 
measure  of  how  the  average  energy  of  the  precipitating 
electrons  changed  when  low-energy  particles  were  present 
in  significant  amounts.  This,  in  turn,  produced  inaccurate 
estimates  of  composition.  Thus  the  technique  presented 
below  is  the  one  that  both  in  theory  and  practice  produced 
results  with  the  least  scatter  and  best  estimate  of  both  the 
energy  and  composition.  A  future  publication  will  review 
the  various  techniques  in  greater  detail  than  is  presented  by 
Hecht  et  al.  [2006]. 

[10]  However,  the  average  energy  of  the  precipitating 
particles  by  itself  is  also  an  important  quantity  to  be  derived 
for  this  study.  This  quantity  can  be  obtained  from  other 
techniques.  Thus  the  average  energy  results  will  be  supple¬ 
mented  by  Doppler  measurements  of  the  OI  (557.7  nm) 
emission  available  at  Poker  Flat,  imaging  riometer  data 
available  at  Poker  Flat,  J  sensor  and  TED  sensor  in  situ  particle 
measurements,  and  TIMED  observations.  The  energy 
deposited  in  an  atmospheric  column  by  particles  and  by  Joule 
heating  will  be  estimated  from  photometer  and  magnetometer 
data.  Because  some  of  these  techniques  are  not  widely  known 
they  will  be  described  in  detail  below. 

2.1.  Ground-Based  Instruments 

2.1.1.  Four  Channel  Photometers  at  Fort  Yukon  and 
Poker  Flat  Alaska 

2. 1.1.1.  Instrument  Description 

[n]  Two  photometers,  each  with  a  four-channel  filter 
wheel,  have  been  deployed  since  2001  at  two  Alaska  sites, 
Poker  Flat  and  Fort  Yukon,  whose  locations  are  given  in 
Table  1.  The  filters  in  each  photometer  are  centered  around 
the  wavelength  for  the  following  auroral  emissions:  the 
(0,  1)  band  of  the  N21N  system  at  427.8  nm  referred  to  as 
blue,  OI  (630.0  nm)  referred  to  as  red,  OI  (844.6  nm),  and 
the  (2,  1)  band  of  the  N2  IP  system  at  871.0  nm.  All  four 
filters  have  bandwidths  of  between  1  and  2  nm  full  width  at 
half  maximum  and  are  typically  sampled  every  10  s.  The 
integration  time  on  each  filter  is  about  1  s.  The  photometers 
look  up  the  local  magnetic  zenith  and  the  instrument  runs 
automatically  from  dusk  to  dawn.  A  version  of  this  system 
has  proven  reliable  and  has  been  used  for  earlier  studies 
[Christensen  et  al .,  1997]. 


[12]  All  emissions  are  corrected  for  backgrounds  by 
subtracting  emissions  measured  during  periods  of  no  aurora. 
Since  auroral  emissions  can  occur  continuously  at  times, 
fixed  backgrounds  of  less  than  100  R  are  subtracted  in 
those  cases.  These  backgrounds  are  derived  from  previous 
measurements.  To  make  sure  that  background  variations  do 
not  add  significantly  to  the  uncertainty  in  the  derived 
results,  most  of  the  analysis  is  restricted  to  periods  of  bright 
aurora  where  the  brightness  of  the  blue  emission  is  above 
250  R. 

2.1. 1.2.  Analysis  Technique 

[13]  The  analysis  technique  to  determine  composition  and 
energetics  has  been  described  in  detail  in  several  previous 
publications  [Hecht  et  al .,  1995;  Christensen  et  al .,  1997; 
Hecht  et  al .,  1999,  2000,  2006]  and  is  based  on  producing  a 
lookup  table  of  predicted  red,  blue,  and  OI  (844.6  nm) 
emissions  using  the  Strickland  electron  transport 
code  [Strickland  et  al. ,  1989].  This  is  done  by  inputting  a 
Gaussian  electron  energy  distribution  of  average  energy  E 
into  the  top  of  a  model  atmosphere  whose  atomic  oxygen 
density  [O]  profile  is  scaled  by  a  parameter  fG.  A  value  of 
one  for  f0  means  that  the  [O]  values  are  those  given  by  the 
model  atmosphere  used  by  the  Strickland  code.  Values  less 
than  one  mean  reduced  [O].  Among  the  many  outputs  of  the 
Strickland  code  are  the  brightnesses  in  Rayleighs  (R)  of  all 
the  auroral  optical  emissions  as  seen  from  the  ground 
discussed  in  this  paper.  Four  different  fG  values  and  seven 
different  energies  form  the  basis  of  our  lookup  table.  From 
the  table  a  plot  of  red/blue  versus  OI  (844.6  nm)/blue  can  be 
used  to  immediately  determine  values  of  E  and  f0  from 
observations  of  these  ratios.  This  latter  parameter,  however, 
only  makes  sense  if  the  details  of  the  model  atmosphere  are 
known.  In  this  work  we  compare  not  only  variations  during 
a  day  but  also  day-to-day  variations  and  the  reader  may  not 
necessarily  have  access  to  how  [O]  varies  in  each  model 
atmosphere  associated  with  every  plot  of  f0.  Furthermore, 
we  also  are  making  comparisons  with  satellite  dayglow 
results  which  report  the  quantity  0/N2  instead  of  f0.  Next, 
we  show  how  the  two  approaches  are  equivalent  and  the 
advantages  of  using  0/N2,  which  we  make  use  of  in  the 
remainder  of  the  paper. 

[14]  To  get  around  this  limitation,  we  introduce  0/N2,  the 
ratio  of  column  densities  referenced  to  a  fixed  N2  column 
density  (similar  to  referencing  a  quantity  to  a  fixed  pressure 
level).  Strickland  et  al.  [1995]  introduced  this  ratio  for 
remote  sensing  applications  with  the  intended  use  of  OI 
(135.6  nm)  and  N2  Lyman-Birge-Hopfield  (LBH)  dayglow 
observations.  The  choice  of  the  reference  column  density  of 
N2  (abbreviated  to  N2  depth)  is  not  critical  but  once 
selected,  must  be  retained  for  convenient  comparisons  of 
derived  0/N2  from  data  set  to  data  set.  The  recommended 


3  of  28 


A01310 


HECHT  ET  AL.:  THERMOSPHERIC  COMPOSITION  ON  17-21  NOV  2003 


A01310 


value  in  the  above  paper  is  1017  cm-2  due  to  slightly  less 
scatter  between  01  (135.6  nm)/LBH  and  [0]/[N2]  than  at 
larger  or  smaller  reference  values  (scatter  plots  were  exam¬ 
ined  from  calculations  utilizing  an  extensive  set  of  model 
atmospheres).  This  recommended  value  typically  occurs 
when  the  integration  as  a  function  of  altitude  which  starts 
at  the  top  of  the  atmosphere  stops  in  the  vicinity  of  135  km. 
Here,  we  also  use  10l7  cm~2  but  emphasize  that  essentially 
the  same  results  follow  from  a  value  deeper  in  the  atmo¬ 
sphere  such  as  1018  cm-2. 

[is]  To  clarify  this,  consider  producing  two  lookup  tables 
from  the  above  table,  one  that  maps  the  fG  grid  to  0/N2 
values  referenced  to,  say,  an  N2  depth  of  1017  cm-2  and  the 
other  to,  say,  1018  cm-  .  The  tables  are  identical  except  for 
their  0/N2  grids.  Use  of  the  tables  with  optical  data  at  a 
given  point  in  time  lead  to  two  different  0/N2  values,  each 
referenced  to  one  of  these  two  depths.  For  any  given 
atmosphere,  we  now  wish  to  know  the  fG  values 
corresponding  to  the  just  derived  optical-based  0/N2  values. 
We  first  perform  the  integrations  to  obtain  the  0/N2  values 
(call  these  model  values)  at  the  two  depths  (these  corre¬ 
spond  to  f0  =  1  for  that  atmosphere).  For  each  of  the 
reference  depths,  the  optical-based  f0  is  just  the  ratio  of 
the  optical-based  0/N2  to  the  model  value.  For  the 
model  atmosphere  used  to  generate  the  tables,  the  two 
derived  0/N2  values,  of  course,  lead  back  to  exactly  the 
same  iQ  value.  For  other  atmospheres,  the  two  f0  values  will 
differ  and  the  question  is  by  how  much.  The  analysis  by 
Strickland  et  al.  [1995]  done  for  the  dayglow  suggests  the 
spread  in  this  difference  is  insignificant  within  the  error 
budget  for  data  and  model  uncertainties  and  is  also  small 
compared  to  the  atmospheric  variability  during  disturbed 
periods  to  be  illustrated  shortly.  The  uncertainty  is  due  to 
the  lack  of  knowledge  as  to  the  correct  model  atmosphere. 
For  column  densities  of  1017  cm-2  they  find  uncertainties  to 
be  around  two  percent  which  is  still  only  around  5  percent 
when  the  column  density  is  1018  cm-2.  These  low  values, 
which  do  not  vary  much  with  column  density,  are  due  to  the 
physically  constrained  smooth  profiles  found  in  model 
atmospheres.  For  an  auroral  case  this  could  be  different. 
However,  we  will  show  that  the  choice  of  a  model  atmo¬ 
sphere  does  not  significantly  affect  our  retrievals. 

[i6]  Note  the  following  about  this  approach:  (1)  All 
auroral  emissions  are  still  calculated  through  the  entire 
model  atmosphere  as  seen  from  the  ground.  (2)  The  only 
difference  from  the  previous  approach  is  that  now  for  every 
model  atmosphere  used  to  produce  a  new  lookup  table,  the 
value  of  i0  of  1  is  replaced  with  the  calculated  0/N2  ratio.  An 
fG  of  0.5  (for  the  same  atmosphere)  will  produce  an  0/N2 
ratio  that  is  half  of  the  value  for  f0  of  1.  (3)  By  using  this 
technique,  however,  one  would  expect  that  the  0/N2 
corresponding  to  an  f0  of  1  for  a  model  atmosphere  during 
a  geomagnetically  quiet  period  would  be  different  than  the 
0/N2  corresponding  to  an  fG  of  1  for  a  geomagnetically 
disturbed  period.  This  will  be  illustrated  in  a  later  section. 
(4)  Finally,  on  the  basis  of  the  the  error  analysis  by 
Strickland  et  al.  [1995]  and  the  above  discussion,  one  would 
expect  to  derive  for  a  given  set  of  photometric  data 
essentially  the  same  0/N2,  independent  of  the  model  atmo¬ 
sphere  used  to  derive  the  lookup  table.  This  will  be  shown 
to  be  the  case  later  in  this  paper. 


[17]  Two  different  model  atmospheres  are  used  here  in  the 
analysis,  both  based  on  the  Naval  Research  Laboratory 
Mass  Spectrometer  and  Incoherent  Scatter  (NRLMSIS) 
model  [Picone  et  al.,  2002]  and  both  run  for  the  location 
of  Poker  Flat.  The  MSIS  series  of  models  are  semiempirical 
models  that  rely  on  physically  constrained  vertical  profiles 
based  on  mass  spectrometer  and  inchorent  scatter  radar  data. 
NRLMSIS  is  an  updated  version  of  MSIS90  [ Hedin ,  1991]. 
For  analysis  on  20  November  and  21  November  the 
NRLMSIS  model  atmosphere  used  is  the  one  appropriate 
to  the  superstorm  conditions  on  20  November.  Thus  the 
FI 0.7  and  F10.7A  (the  81  d  average)  are  taken  as  155.1  and 
138.3,  respectively,  and  the  3  h  a p  mode  is  used  appropriate 
to  a  universal  time  (UT)  of  15  h.  The  ap  values  are  150.00 
(daily),  300.00  (0-3  h),  179.00(3-6  h),  94.00  (6-9  h), 
94.00  (9-12  h),  12.50  (12-33  h),  and  19.38  (33-59  h), 
where  the  first  is  the  daily  value  and  the  others  are  the  sum 
of  the  a p  values  over  the  period  prior  to  1 500  UT.  For  this 
model  atmosphere  the  column  0/N2  is  0.449.  Although  this 
time  period  is  during  a  superstorm,  for  ease  of  referral  this  is 
called  the  “disturbed”  atmosphere. 

[is]  For  analysis  on  17  November  and  19  November  the 
NRLMSIS  model  atmosphere  used  is  the  one  appropriate  to 
the  less  disturbed  conditions  found  around  1000  UT  on 
19  November.  Thus  the  FI 0.7  and  F10.7A  are  taken  as 
144.3  and  138.3,  respectively,  and  again  the  3  h  ap  mode  is 
used  appropriate  to  a  universal  time  of  15  h.  The  a p  values 
are  1 1.50  (daily),  9.00  (0-3  h),  15.00(3-6  h),  9.00  (6-9  h), 
9.00  (9-12  h),  25.88  (12-33  h),  and  39.00  (33-59  h), 
where  the  first  is  the  daily  value  and  the  others  are  the  sum 
of  the  a p  values  over  the  period  prior  to  1000  UT.  For  this 
model  atmosphere  the  column  0/N2  is  0.91.  This  is  referred 
to  as  the  “quiet”  atmosphere. 

[19]  Figures  la  and  lb  show  the  library  of  the  Strickland 
model  runs  of  red/blue  versus  OI  (844.6  nm)/blue  for  these 
two  atmospheres  based  on  Gaussian  energy  distributions 
with  assumed  low-  and  high-energy  tails  [Strickland  et  al, 
1993].  This  will  be  referred  to  as  an  enhanced  Gaussian.  For 
a  Gaussian  without  tails  the  characteristic  energy  Eq  is  equal 
to  the  energy  of  the  distribution  at  the  peak  and  to  the 
average  energy  E  of  the  distribution.  For  the  particle 
distributions  used  here  with  tails  the  characteristic  and 
average  energies  are  the  same  to  within  10%  and  thus  only 
E  will  be  used  in  the  remainder  of  the  paper.  Similar  plots 
have  been  presented  previously  [e.g.,  see  Hecht  et  al, 
1989],  which  allow  an  easy  visualization  of  how  0/N2 
and  E  depend  on  these  ratios.  Solid  lines  show  regions  of 
constant  0/N2  and  dotted  lines  show  regions  of  constant  E. 

[20]  This  analysis  of  the  photometric  data  assumes  that 
the  incoming  electron  distribution  is  characterized  by  the 
enhanced  Gaussian.  Gaussians  appear  to  be  a  reasonable 
choice  based  on  prior  observations  during  large  geomag¬ 
netic  storms  and  during  moderate  substorm  activity  [Hecht 
et  al,  1991,  1999].  While  Maxwellians  sometimes  are 
found,  these  tend  to  occur  during  diffuse  aurora  which 
usually  have  small  energy  fluxes.  As  will  be  shown  later, 
there  is  some  uncertainty  involved  in  this  assumption.  On 
average  for  the  data  considered  here  Maxwellians  increase 
the  derived  0/N2  by  25%  and  the  derived  E  by  10%.  The 
reason  for  this  increase,  especially  in  0/N2  is  partly 
explained  as  follows.  Strickland  model  results  show  that 
the  centroid  of  emission  for  a  Maxwellian  distribution  is 
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Figure  1.  Plots  of  the  derived  0/N2  (solid  line)  and  E  (dotted  line)  values  using  the  Strickland  model 
and  various  algorithms  (see  text),  showing  (a)  red/blue  versus  OI  (844.6  nm)/blue  and  “quiet” 
atmosphere,  (b)  red/blue  versus  OI  (844.6  nm)/blue  and  “disturbed”  atmosphere,  and  (c)  red/blue 
versus  OI  (557.7  nm)/blue  and  “disturbed”  atmosphere.  Note  that  0/N2  of  1.35  is  not  labeled.  Also 
shown  is  (d)  the  derived  E  values  as  a  function  of  OI  (557.7  nm)  temperature  with  solid  line  denoting 
“quiet”  atmosphere  and  dotted  line  denoting  “disturbed”  atmosphere. 


lower  in  altitude  (where  both  distributions  are  constrained 
to  replicate  the  observations).  In  order  to  replicate  the 
observed  01(844. 6nm)/blue  ratio  in  this  situation,  the 
Gaussian-based  O  density  profile  must  be  increased  when 
switching  to  a  Maxwellian  distribution  to  compensate 
for  the  reduction  in  this  ratio  by  moving  the  emission 
centriod  downward  in  the  atmosphere.  However,  there  is 
also  a  technique  for  determining  whether  Gaussians  are 
appropriate  [ Christensen  et  al ,  1987;  Strickland  et  al. , 
1994]  and  we  will  apply  it  at  the  end  of  the  paper  to  show 
that  the  Gaussian  assumption  is  reasonable. 

[21]  Now  consider  what  is  inferred  from  this  technique  in 
an  extreme  example  when  equal  energy  fluxes  from  a 


10  keV  and  0.5  keV  enhanced  Gaussian  energy  distributions 
are  present  for  an  0/N2  ratio  of  0.23.  The  inferred  E  value 
will  be  keV  while  the  0/N2  will  stay  nearly  the  same. 
The  inferred  E  value  of  around  1  keV  is  technically  correct, 
although  its  use  to  derive,  say,  the  electron  density  profile, 
would  produce  a  significant  error.  This  is  due  to  the 
presence  of  high-energy  particles  going  undetected. 

[22]  Furthering  this  example,  consider  the  sensitivity  of 
the  two  derived  quantities  when  that  distribution  is 
enhanced  or  deficient  in  high-  or  low-energy  “tails.”  Both 
the  red/blue  and  OI  (844.6  nm)/blue  ratios  show  sensitivity 
to  the  presence  of  low-  (below  1  keV)  and  high-  (above 
10  keV)  energy  electrons  as  these  ratios  change  with 
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Table  2.  Uncertainties  in  the  Analysis  of  Ground-Based  Data 


Title 

%  for  0/N2 

%  for  E 

Counting  Statistics 

less  than  5 

less  than  5 

Model  Atmosphere 

5 

10 

Reference  Depth 

5 

5 

Energy  Distribution 

(Maxwellian  or  Gaussian) 

25 

10 

Absolute  Calibration 

20 

20 

Relative  Calibration 

10 

10 

Total  Absolute 

32 

23 

Total  Relative 

28 

16 

Total  Relative  with 
known  Gaussian 

12 

12 

E.  However,  using  photometric  data  to  derive  E  and  0/N2 
will  show  more  sensitivity  to  low-energy  electrons.  For 
example,  consider  where  a  case  of  moderate  blue  emission, 
on  the  order  of  1000  R,  and  both  ratios  are  at  least  0.3  or 
larger.  The  E  value  will  be  less  than  5  keV  and  both  the  red 
and  OI  (844.6  nm)  emissions  will  be  significant  (300  R  or 
much  greater).  In  this  case  there  will  be  a  high  sensitivity  to 
low-energy  electrons.  In  contrast,  if  both  ratios  are  around 
0. 1  or  less,  E  will  be  above  10  keV.  This  means  both  the  red 
and  OI  (844.6  nm)  emission  will  be  small  (on  the  order  of 
100  R  or  below).  Thus  the  low  ratios  will  be  more  affected 
by  measurement  uncertainties  and  background  corrections 
in  the  red  and  OI  (844.6  nm)  emissions.  In  addition, 
independent  of  the  blue  emission  brightness,  this  greater 
sensitivity  to  low-energy  electrons  is  reflected  in  the  com¬ 
position  measurement  as  the  separation  of  the  lines  of 
constant  0/N2  increase  toward  lower  energy.  Thus  the 
increasing  sensitivity  of  both  ratios  to  low-energy  electrons 
means  that  a  small  enhancement  of  the  low-energy  compo¬ 
nent  can  be  measured  more  easily  than  a  high-energy 
enhancement  resulting  in  a  tendency  to  obscure  the  presence 
of  high-energy  particles  while  emphasizing  the  presence  of 
a  low-energy  component.  Nevertheless,  for  all  but  very  low 
electron  fluxes,  the  E  values  and  the  0/N2  results  are  not 
much  affected  by  uncertainties  in  the  shape  of  the  assumed 
Gaussian  energy  distribution. 

[23]  The  electron  flux  Q  can  be  derived  from  the  blue 
intensity.  Nominally,  1  mW/m2  (1  erg-cm“2-s_1)  of  elec¬ 
trons  produces  about  250  R  of  blue.  Since  the  proton  flux  is 
usually  a  small  fraction  of  the  electron  flux  and  the  emission 
efficiency  for  protons  and  electrons  is  similar  [Strickland  et 
al.,  1993],  any  difference  between  the  species  is  ignored. 
The  derived  Q  reported  in  this  work,  called  Q(particle)  is, 
except  where  noted,  the  total  electron  and  proton  flux.  All 
plots  are  interpolated  onto  1  min  grids. 

[24]  In  previous  studies  we  have  estimated  the  uncertainty 
to  be  typically  less  than  30%  and  often  better  than  10% 
[Hecht  et  al .,  2006].  However,  with  the  comparisons  we 
show  later  we  can  better  quantify  this.  Table  2  lists  the 
major  sources  of  error  and  the  associated  effects  on  the 
derived  0/N2  and  E.  For  absolute  uncertainties,  all  must  be 
considered,  while  for  measurements  of  relative  changes  the 
absolute  calibration  is  replaced  by  the  relative  calibration. 
Note,  though,  that  the  total  uncertainties  are  actually  for  a 
single  point  as  the  calibration  uncertainties  apply  to  all  data 
and  do  not  vary  from  measurement  to  measurement.  Thus 
the  quoted  total  uncertainties  represent  an  upper  limit.  Since 
a  knowledge  of  the  energy  distribution  is  the  major 


unknown  it  appears  that  in  a  relative  sense  we  can  obtain 
0/N2  to  about  30%  and  E  to  about  15%  if  we  have 
no  knowledge  of  whether  a  Gaussian  is  a  reasonable 
representation.  If  such  knowledge  is  available,  then  these 
uncertainties  decrease  to  less  than  15%. 

[25]  However,  during  dynamic  aurora,  the  error  in  a 
derived  0/N2  or  E  could  exceed  these  numbers.  This  occurs 
because  during  the  5  s  observation  period  the  aurora 
changes  so  that  the  filter  ratios  are  no  longer  a  valid 
representation  to  be  used  in  this  analysis.  Owing  to  this 
uncertainty,  we  will  limit  interpretation  in  the  present  to 
long  period  changes  that  occur  over  tens  of  minutes  rather 
than  short-term  fluctuations  over  a  few  minutes  or  so. 
2.1.2.  Meridian  Scanning  Photometers  at  Poker  Flat 
Alaska,  Gillam  and  Pinnawa  Canada 

2.1.2.1.  Instrument  Description 

[26]  The  Poker  Flat  Observatory  [Lummerzheim  et  al ., 
1990]  as  well  as  two  observatories  in  central  Canada 
at  Pinnawa  and  Gillam  (see  Table  1)  have  Meridian  Scan¬ 
ning  Photometers  (MSP)  that  measure  OI  (557.7  nm),  red 
and  H-beta  (486.1  nm).  At  Poker  the  MSP  measures  blue 
while  the  MSPs  at  Pinnawa  and  Gillam  measure  the  (0,2) 
band  of  the  N21N  system  at  470.9  nm.  The  blue  and  the 
470.9  nm  bands  are  within  the  same  system  and  their  ratios 
are  constant  [Vallance  Jones ,  1974].  These  measurements 
are  made  along  the  magnetic  meridian  with  data  samples  up 
the  magnetic  zenith  occurring  approximately  once  a  minute. 
For  the  Poker  MSP,  the  data  have  not  been  recently 
calibrated.  Thus  in  the  analysis  below  the  Poker  MSP  data 
will  only  provide  information  on  H-beta.  The  H-beat  data, 
even  in  a  relative  sense,  provide  insight  as  to  correlations 
between  the  occurrence  of  proton  aurora  and  changes  E 
obtained  using  the  photometer  data  described  above. 
Furthermore,  since  the  derived  0/N2  is  sensitive  to  the 
absolute  calibration,  only  the  E  results  will  be  presented 
for  the  Canadian  stations. 

2.1.2.2.  Analysis  Technique 

[27]  The  H-beta  measurements  allow  a  measure  of  proton 
energy  flux  since  ~50  to  100  R  of  H-beta  is  about  1  mW/m2 
[Vallance  Jones ,  1974].  The  range  exists  because  the  flux 
depends  on  the  average  energy  of  the  protons  [Vallance 
Jones ,  1974;  Strickland  et  al ,  1993].  Recent  studies  of 
satellite  particle  data  [Coumans  et  al .,  2002]  suggest 
average  energy  values  well  above  20  keV  for  protons,  and 
thus  we  adopt  a  nominal  H-beta  yield  of  50  R  for  H-beta 
emission  for  1  mW/m2  of  incoming  protons. 

[28]  Although  the  emission  efficiency  of  OI  (557.7  nm) 
may  be  higher  for  proton  excitation  than  for  electron 
excitation,  this  remains  uncertain.  Since  previous  workers 
[Vallance  Jones ,  1974]  have  found  similar  efficiencies  for 
protons  and  electrons,  we  will  adopt  (as  we  did  for  the  N2 
and  other  OI  emissions)  the  same  efficiency  for  proton  and 
electron  excitation. 

[29]  The  Strickland  model  includes  the  auroral  excitation 
and  quenching  of  O  ('S)  responsible  for  OI  (557.7  nm)  and 
as  discussed  by  Strickland  et  al  [2000],  the  two  main 
sources  are  direct  electron  impact  excitation  of  O,  and  the 
quenching  of  the  N2(A)  state  by  O.  The  latter  mechanism 
dominates  in  the  lower  thermosphere  where  the  bulk  of 
excitation  occurs  under  active  conditions.  At  low  0/N2 
values  <0.5,  OI  (557.7  nm)  increases  with  increasing  0/N2. 
The  quenching  of  N2(A)  in  this  situation  by  O  has  a  weak 
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impact  on  the  N2(A)  density.  However,  at  the  larger  E 
values  (e.g.,  >5  keV)  where  most  of  the  excitation  is  in 
the  lower  thermosphere,  01  (557.7  nm)  becomes  insensitive 
to  increases  in  0/N2  as  it  extends  beyond  ~1 .  This  satura¬ 
tion  effect  arises  from  an  increase  in  the  O  density  to  the 
point  of  reducing  the  density  of  N2(A)  through  quenching. 
Thus  the  volume  excitation  rate  of  O  (!S)  for  this  reaction 
(proportional  to  the  product  of  these  densities)  becomes 
effectively  constant. 

[30]  Figure  lc  shows  a  plot  of  01  (557.7  nm)/blue  versus 
red/blue  as  functions  of  E  and  0/N2  for  the  “disturbed” 
atmospheres  that  is  used  for  analyzing  the  Canadian  MSP 
data  on  20  November.  Note  that  the  derivation  of  E  depends 
greatly  on  the  red/blue  ratio.  Even  if  an  uncertain  calibration 
produces  inaccuracies  in  the  absolute  value  of  E,  relative 
changes  in  E,  which  are  the  intended  use  of  these  data, 
remain  accurate. 

2.1.3.  Scanning  Doppler  Imaging  Fabry-Perot  at 
Poker  Flat,  Alaska 

2. 1.3.1.  Instrument  Description 

[31]  The  Scanning  Doppler  Imaging  Fabry-Perot  devel¬ 
oped  at  Poker  Flat  Observatory,  Alaska  [Conde  and  Smith , 
1998]  measures  the  line  profile  of  01  (557.7  nm)  allowing 
an  inference  of  E  from  the  measured  temperature.  The 
instrument  actually  produces  a  map  of  the  temperature  over 
the  whole  sky.  For  this  study  the  pixels  closest  to  the 
magnetic  zenith  are  used. 

2. 1.3. 2.  Analysis  Technique 

[32]  As  noted  by  Hecht  et  al.  [2006],  this  technique 
provides  an  alternative  method  for  deriving  E  and  0/N2. 
The  Strickland  model  includes  all  relevant  cross  sections, 
chemical  reactions,  and  quenching  that  produce  01 
(557.7  nm)  and  thus  allows  the  01  (557.7  nm)  Doppler 
temperature  to  be  predicted  for  a  given  E. 

[33]  Figure  Id  shows  the  relationship  between  measured 
temperature  and  E  for  the  geomagnetically  “quiet”  and 
“disturbed”  atmospheres.  As  this  technique  depends  on 
the  OI  (557.7  nm)  temperature,  it  is  not  very  sensitive  to 
low-energy  electrons  since  most  of  the  OI  (557.7  nm) 
emission  occurs  in  the  lower  E  region.  For  the  hypothetical 
example  above  of  equal  0.5  and  10  keV  electron  energy 
fluxes,  the  inferred  E  value  would  be  ~7  keV  assuming  the 
“quiet”  atmosphere  shown  in  Figure  Id.  The  uncertainties 
in  the  zenith  temperatures  are  on  the  order  of  1%,  which 
translates  into  less  than  0.1  keV  error. 

[34]  There  are  potential  problems  at  energies  above 
10  keV.  First,  at  very  high  energies  well  above  10  keV, 
the  emission  comes  from  an  altitude  below  100  km  where 
OI  (557.7  nm)  begins  to  be  severely  quenched.  This  is 
reflected  in  the  results  for  the  “quiet”  atmosphere  in  which 
the  Doppler  temperature  shows  only  small  temperature 
changes  as  E  approaches  25  keV.  Second,  dynamical  effects 
such  as  gravity  waves  or  tides  are  not  included  in  the  model 
atmosphere  and  these  can  distort  the  temperature  profile  and 
therefore  the  inferred  temperature.  This  is  especially  signif¬ 
icant  when  high-energy  electrons  are  present  that  penetrate 
to  around  100  km,  a  region  known  for  strong  wave  activity 
and  where  significant  OI  (557.7  nm)  emission  can  occur.  In 
addition  to  these  two  problems,  there  is  the  possibility  that 
the  model  atmosphere  may  simply  not  represent  the  actual 
atmospheric  temperature  profile  accurately.  An  example  of 
this  possibility  is  given  in  a  later  section. 


2.1.4.  Magnetometers  at  Poker  Flat  and 
Fort  Yukon,  Alaska 

2.1.4.1.  Instrument  Description  and  Analysis  Technique 

[35]  Following  Duboin  and  Kamide  [1984]  and  ignoring 
the  effect  of  neutral  winds  [Thayer  et  al ,  1995;  Thayer , 
1998],  the  amount  of  electromagnetic  energy  deposited  in 
a  volume  per  unit  time  is  the  Joule  heating  rate,  q,,  in 
W/m3.  The  height-integrated  Joule  heating  rate,  Q j  in  units 
of  mW/m2,  equals  the  rate  at  which  electromagnetic  energy 
is  deposited  into  an  atmospheric  column  [Duboin  and 
Kamide ,  1984;  Thayer  et  al,  1995].  Ignoring  neutral  winds, 
Qj  calculated  following  Duboin  and  Kamide  [1984]  is 
equivalent  in  units  to  a  flux  [Thayer  et  al ,  1995]  and 
represents  the  electromagnetic  energy  deposited  into  the 
ionosphere  as  thermal  energy  [Thayer  et  al ,  1995].  How¬ 
ever,  neutral  winds  have  an  effect  on  this  quantity  as 
discussed  by  Thayer  et  al  [1995]  and  Thayer  [1998]. 
However,  even  in  those  studies  for  the  most  part,  Qj 
calculated  ignoring  winds  often  provided  a  good  proxy  for 
the  height-integrated  Joule  heating  rate  including  neutral 
wind  effects. 

[36]  Duboin  and  Kamide  [1984]  suggest  that  Qj  can  be 
estimated  from  geomagnetic  field  data  obtained  from 
ground-based  magnetometers.  Magnetometers  are  present 
at  both  Alaska  sites  and  these  data  have  been  used  in 
previous  studies  to  estimate  Qj  [Duboin  and  Kamide , 
1984;  Christensen  et  al ,  1997;  Hecht  et  al,  1995,  2000, 
2006].  Often  Q j  well  exceeds  Q  (particle)  and  thus  Qj  may 
be  the  primary  cause  of  vertical  winds  and  subsequent 
composition  change  [e.g.,  see  Christensen  et  al,  1997]. 

[37]  To  obtain  Qj,  we  use  the  AH  component  from  a 
magnetometer  (i.e.,  the  change  of  the  H  component  from 
nominally  nonauroral  conditions  [Duboin  and  Kamide, 
1984]).  The  value  of  (AH)2  in  units  of  (nT)2  can  be  used 
when  multiplied  by  a  constant,  8  x  10“5  mW/((m2XnT2)), 
to  estimate  Q j,  although  this  approximation  may  be 
inaccurate  on  occasion.  In  particular  from  their  study,  which 
took  place  in  April  1978,  Duboin  and  Kamide  [1984]  stated 
that  the  above  constant  was  appropriate  during  nighttime 
periods  (after  about  0800  UT  at  Poker  Flat)  when  AH  was 
negative,  indicative  of  a  westward  electrojet.  During  the 
early  evening  prior  to  0800  UT  when  AH  was  positive  and 
indicative  of  an  eastward  electrojet,  the  constant  was  found 
to  be  over  10  times  larger  at  105  x  10~5  mW/((m2XnT2)). 
For  the  vast  majority  of  data  presented  in  this  study,  the 
lower  value  is  used,  although  we  have  used  the  higher 
constant  where  appropriate. 

[38]  In  a  previous  study,  comparisons  have  been  made 
during  nighttime  auroral  activity  between  Q j  calculated 
from  radar  data  at  Kangerlussuaq,  Greenland  and  the 
approximation  to  Qj  used  in  this  paper  [Hecht  et  al, 
2000].  The  use  of  the  lower  constant  approximation  for 
periods  of  negative  AH  was  found  to  provide  a  good 
approximation  to  Qj  independent  of  the  original  analysis 
by  Duboin  and  Kamide  [1984]. 

2.1.5.  Imaging  Riometer  at  Poker  Flat,  Alaska 

2.1.5.1.  Instrument  Description 

[39]  A  riometer  observes  ionospheric  absorption  of  radio 
noise  (called  cosmic  noise  absorption;  CNA),  which  is  used 
to  infer  the  lower  ionospheric  electron  density  variation.  A 
riometer  generally  operates  at  frequencies  of  20  to  50  MHz. 
Observed  from  the  ground,  cosmic  radio  noise  is  partially 


7  of  28 


A01310 


HECHT  ET  AL.:  THERMOSPHERIC  COMPOSITION  ON  17-21  NOV  2003 


A01310 


absorbed  by  the  ionosphere  in  the  80-90  km  altitude  range 
and  the  amount  of  the  radio  absorption  will  vary  depending 
on  ionospheric  conditions.  Disturbances  in  electron  density 
in  the  lower  ionosphere  can  be  measured  by  observing  the 
variation  of  this  cosmic  noise  intensity  on  the  ground.  In 
polar  regions  aurora  is  one  of  the  events  which  can  cause 
CNA  variation. 

[40]  The  imaging  riometer  at  Poker  Flat  has  been  in 
continuous  operation  since  October  1995  and  it  consists 
of  a  1 6  x  1 6-element  square  array  antenna,  a  Butler  matrix, 
16  receivers,  and  personal  computers  for  system  control. 
The  imaging  riometer  scans  an  antenna  beam  in  about  200 
directions  within  1  s  to  monitor  cosmic  radio  noise  variation 
at  38.2  MHz  over  a  wide  field  of  view  of  the  radio  sky.  It 
observes  lower  ionospheric  disturbances  over  a  400-km  by 
400-km  field  of  view  at  the  90-km  height  with  a  spatial 
resolution  of  about  7°  (11  km  around  the  zenith),  a  sensi¬ 
tivity  better  than  0.1  dB,  and  a  time  resolution  of  1  s 
[Murayama  et  al,  1997;  Mori  et  al. ,  2004]. 

2.1.5.2.  Analysis  Technique 

[41]  For  data  analysis  in  this  work,  we  use  the  mean  value 
of  CNA  for  the  two  beams  closest  to  the  geomagnetic 
zenith.  The  presence  of  CNA  suggests  that  high-energy 
electrons  are  involved,  since  electron  energies  greater  than 
25  keV  are  required  to  penetrate  down  to  altitudes  where 
CNA  occurs.  Kosch  et  al  [2001]  and  Mori  et  al  [2004] 
have  recently  discussed  the  utility  of  using  imaging 
riometers  to  determine  E  based  on  the  relationship  that  E 
is  proportional  to  Ab/Q(particle)° 5  where  Ab  is  the  mea¬ 
sured  riometer  absorption  in  decibels.  Thus  combining 
Q(particle)  values  derived  from  the  photometer  data  and 
Ab  from  the  riometer  data  allow  another  measure  of  E, 
albeit  one  sensitive  to  high-energy  but  not  low-energy  tails. 
Because  of  the  proportionality,  an  unknown,  arbitrary  fixed 
scaling  factor  must  be  applied  and  this  is  done  by  scaling 
the  riometer  results  by  a  factor  of  20.  This  factor  was  chosen 
to  provide  a  reasonable  agreement  with  the  SDI-derived  E 
values  on  17  November.  The  use  of  a  scaling  factor, 
however,  only  affects  the  absolute  comparisons  and  not 
the  relative  comparisons  with  other  instruments. 

[42]  Continuing  the  same  example  used  for  the  previous 
techniques,  the  riometer  would  infer  an  £  of  7  keV  as  Q 
would  double  and  Ab  would  essentially  stay  the  same.  We 
note  that  the  radio  star  (Cassiopeia-A)  is  located  just  on 
the  geomagnetic  zenith  (beam  number  [N9E8])  at  0600- 
0700  UT  during  17-21  November  2003.  This  means  CNA 
around  0500-0800  UT  may  be  affected  by  the  scintillations 
and  the  large  source  signal  of  this  radio  star. 

2.2.  Space-Based  Instruments 

2.2.1.  POES  Particle  Detectors 

2.2.1. 1.  Instrument  Description  and  Analysis  Technique 

[43]  The  NOAA-16  and  NOAA-17  POES  satellites  each 
include  a  Total  Energy  Detector  (TED)  instrument  that 
measures  the  proton  and  electron  flux  in  a  number  of 
channels  from  50  eV  to  20  keV.  Unfortunately,  not  all  data 
are  sent  down  and  the  data  used  in  this  study  are  the  ratios 
of  the  omnidirectional  electron  energy  flux  below  1  keV  to 
the  total  electron  energy  flux  and  the  ratio  of  the  proton  to 
electron  energy  flux.  Because  much  of  the  proton  flux 


occurs  above  30  keV  [Coumans  et  al ,  2002],  we  also  add 
to  the  TED  proton  flux  that  measured  by  the  Medium 
Energy  Proton  and  Electron  Detector  (MEPED)  (part 
of  the  MEPED/TED  package),  which  measures  proton 
energies  from  30  keV  to  1000  keV.  In  November  2003 
the  Local  Time  Ascending  Node  (LTAN)  for  NOAA-16  was 
1409  LST  and  for  NOAA-17  was  2219  LST.  LTAN  is  the 
local  solar  time  (LST)  at  the  subsatellite  point  when  the 
satellite  crosses  the  equator  northbound.  Typical  times  for 
passing  over  or  close  to  a  polar  northern  hemisphere  site 
such  as  Poker  Flat  Alaska  at  night  will  be  quite  different  for 
the  different  satellites.  On  20  November  2003  a  NOAA-17 
northbound  pass  came  within  3°  of  Poker  Flat  at  0702  UT 
(2102  LST)  during  the  early  evening.  A  NOAA-16  south¬ 
bound  pass  also  came  within  3°  of  the  same  site  at  1311  UT 
(0327  LST)  in  the  postmidnight  sector. 

2.2.2.  DMSP  Particle  Detectors 

2.2.2. 1.  Instrument  Description  and  Analysis  Technique 

[44]  The  DMSP  F15  and  F16  satellites  each  carry  a 
particle  detector,  J/4  on  FI 5  and  J/5  on  FI 6.  FI 5  has  an 
LTAN  of  about  2110  LST  while  F 1 6  has  an  LTAN  of  about 
1954  LST.  On  20  November  2003  the  closest  FI 5  and 
FI 6  nighttime  passes  to  Poker  Flat  occurred  near  the 
Alaska-Canada  border  at  0522  and  0350  UT,  respectively. 
Both  were  between  sunset  and  midnight  local  time.  The  two 
J  sensors  differ  in  that  the  J/5  sensor  has  a  wider  field-of- 
view  and  can  obtain  direction  information  for  precipitating 
particles,  while  the  J/4  sensor  looks  only  into  the  zenith 
direction.  For  the  purposes  of  this  work,  the  J/5  and  J/4  data 
both  are  presented  as  flux  integrated  over  the  downward 
hemisphere,  assuming  the  flux  is  independent  of  pitch 
angle.  Each  J  sensor  has  19  energy  channels  for  electrons 
and  another  19  energy  channels  for  protons,  covering  the 
30  eV  to  30000  eV  energy  range.  In  this  study  we  sum  these 
flux  values  to  form  three  energy  channels  from  30  to  100  eV, 
100  to  1000  eV,  and  above  1000  eV.  This  allows  a  more 
precise  measurement  of  the  contribution  of  particles  with 
energies  below  100  eV,  which  are  especially  important  in 
the  production  of  the  rare  type  A  red  aurora  [Robinson  et 
al,  1985]. 

2.2.3.  GUVI 

2. 2.3.1.  Instrument  Description 

[45]  One  of  the  four  instruments  aboard  the  NASA 
TIMED  satellite  is  the  Global  Ultraviolet  Imager  (GUVI), 
an  imaging  spectrograph  capable  of  making  global  maps  in 
five  FUV  colors:  Lyman-alpha  at  121.6  nm,  OI  at  130.4  nm 
and  135.6  nm,  and  two  portions  of  the  N2  Lyman-Birge- 
Hopfield  (LBH)  band  system,  LBH  short  wavelength 
(designated  LBHs)  at  141-153  nm,  and  LBH  long  wave¬ 
length  (designated  LBH1)  at  167  to  181  nm.  While  the  pixel 
sizes  are  quite  small  on  the  ground  (less  than  10  x  10  km 
for  nadir  pixels),  in  this  work  they  are  binned  to  approxi¬ 
mately  35  by  35  km  to  increase  the  sensitivity.  A  full 
description  of  the  instrument  is  given  by  Christensen  et 
al  [2003]. 

2.2.3.2.  Analysis  Technique 

[46]  Ratios  of  the  maps  of  the  OI  (135.6  nm)  to  LBHs 
colors  provide  a  measure  of  0/N2  when  these  data  are 
obtained  from  the  dayglow  [Strickland  et  al,  1995,  2004]. 
The  quiet-time  analysis  by  Strickland  et  al  [2004]  has  been 


8  of  28 


A01310 


HECHT  ET  AL.:  THERMOSPHERIC  COMPOSITION  ON  17-21  NOV  2003 


A01310 


16.0  16.5  17.0  17.5  18.0  18.5  19.0  19.5  20.0  20.5  21.0  21.5  22.0 
Time(UT  Day  of  November  2003) 


Figure  2.  Plots  of  various  quantities  as  a  function  of  UT  from  the  period  of  0000  UT  on  16  November 
2003  to  2400  UT  on  21  November  2003.  Poker  Flat  optical  data  are  only  available  during  the  night. 
Shown  are  (a)  Kp  and  (b)  the  fluence,  F j  over  Poker  Flat  derived  by  integrating  Qj  heating  rate  over  the 
preceding  hour.  For  the  6  d  the  daily  electromagnetic  energy  fluence  (derived  by  integrating  Qy  over  the 
entire  24  h)  is  shown  as  the  series  of  numbers  in  the  upper  part  of  the  panel.  Also  shown  is  (c)  Q(particle), 
designated  here  as  Qp ,  over  Poker  Flat  derived  from  the  blue  photometer  data;  (d)  the  ratio  of  the  Qp 
to  Qt  over  Poker  Flat;  (e)  Poker  Flat  Riometer  absorption  (Ab),  in  decibels,  from  the  imaging 
riometer;  and  (f)  the  average  energy  of  the  precipitating  particles  derived  from  the  Poker  Flat 
riometer  and  photometer  data. 


extended  for  this  November  2003  storm  period  by  Meier 
et  al.  [2005]  and  Crowley  et  al  [2006].  However,  a 
variation  of  the  ground-based  photometric  (red/blue  versus 
01(844.6  nm)/blue)  technique  discussed  above  can  also  be 
used  in  the  auroral  zone  on  the  nightside  to  obtain  E  and 
Q(particle).  The  ratio  of  LBHs  to  LBH1  gives  E  while  either 
one  of  the  three  channel  brightness  gives  Q(particle). 
However,  the  E  algorithm  is  not  valid  for  values  below 
2  keV. 

3.  Results 
3.1.  Overview 

[47]  A  superstorm  occurred  from  about  0800  UT  on 
20  November  2003  to  ~0800  UT  on  21  November  2003. 
Although  the  day  before  the  storm  was  geomagnetically 
quiet,  a  few  days  before  that  there  was  some  geomagnetic 


storm  activity.  To  provide  a  context  for  the  period  from 
0000  UT  on  16  November  to  2400  UT  on  21  November, 
Figure  2  shows  the  index,  the  total  energy  deposited  in 
an  atmospheric  column  due  to  Joule  heating  with  Qj 
integrated  over  the  previous  hour  from  the  Poker  Flat 
magnetometer  data,  Q(particle)  from  the  blue  emission, 
the  ratio  of  Q(particle)  to  the  sum  of  Q(particle)  and  Qy, 
the  Poker  Flat  imaging  riometer  absorption  (Ab),  and  E 
derived  from  combining  the  riometer  Ab  and  Q(particle) 
data  as  described  above  with  the  factor  of  20  scaling  factor 
applied  as  discussed  earlier. 

3.1.1.  Kp 

[48]  The  K p  data  in  Figure  2a  shows  that  significant 
magnetic  activity  occurred  on  16  November  and  17  November, 
which  decreased  on  18  November  and  19  November.  The 
lowest  value  of  Kp  occurred  around  2400  UT  on 
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19  November.  A  large  geomagnetic  storm  occurred  on 

20  November  with  Kp  values  of  8.7.  The  storm  activity 
was  still  present  through  the  first  third  of  21  November  with 
activity  decreasing  throughout  the  rest  of  that  day. 

3.1.2.  Electromagnetic  Energy  Fluence 

[49]  For  a  typical  substorm  time  interval  of  1  h,  the  total 
amount  of  electromagnetic  energy  deposited  in  a  column 
over  that  time  interval  is  derived  by  integrating  Qy  over  the 
preceding  1  h.  Since  Qy  has  the  same  units  as  a  flux  and  the 
time  integral  of  a  flux  is  a  fluence,  in  this  paper  this  integral 
will  be  referred  to  as  the  fluence,  F y.  The  advantage  of 
plotting  F j  rather  than  Qy  is  that  geomagnetic  storm  events 
of  equal  Qy  will  differ  in  Fy  if  one  event  lasts  longer  than  the 
other.  A  plot  of  Fy  also  shows  the  potential  relation  between 
composition  change  and  energy  deposition  in  the  column. 
This  allows  comparison  between  the  energy  deposited  in  a 
column  by  Joule  heating  and  particle  heating. 

[50]  In  evaluating  plots  of  F y,  the  following  should 
be  noted:  if  the  substorm  or  storm  event  lasts  less  than  1  h 
(assuming  a  constant  Qy )  then  the  total  time  when  Fy  is 
above  zero  will  be  1  h  or  greater.  If  the  event  lasts  1  h  then 
the  peak  in  Fy  occurs  at  the  end  of  the  event.  For  events 
lasting  less  than  1  h,  the  highest  value  in  Fy  represents  the 
column  energy  deposition  during  the  event.  If  the  event  lasts 
longer  than  1  h,  the  peak  in  Fy  represents  a  lower  limit  to  the 
energy  deposition. 

[51]  The  data  in  Figure  2b  are  from  the  Poker  Flat 
magnetometer  and  as  expected  typically  show  a  peak  in 
the  postmidnight  sector.  The  peak  values  are  larger  on  lb- 
17  November  than  on  18- 19  November.  The  largest  values 
occur  on  20  November  during  the  storm.  Interestingly,  low 
values  occur  on  21  November. 

[52]  Also  shown  as  numbers  are  the  total  electromagnetic 
energy  deposited  in  a  column  for  each  day  (F Jday),  derived  by 
integrating  Qy  over  the  entire  24  h  period.  These  results  allow 
a  comparison  of  energy  deposition  due  to  the  longer-lasting 
storm  events.  These  values  scale  with  Kp  except  as  noted  on 

2 1  November  where  F Jday  is  the  lowest  of  this  period. 

3.1.3.  Q(particle) 

[53]  The  data  in  Figure  2c  are  derived  from  the  Poker  Flat 
blue  emission  available  only  at  night.  Q(particle)  values  are 
larger  on  16-17  November  than  on  18-19  November,  as 
expected.  While  the  Q(particle)  values  on  20  November 
appear  comparable  to  those  of  the  other  storm  nights  lb- 
17  November,  Q(particle)  on  21  November  is  low.  Note  that 
for  comparisons  with  F y  shown  later,  we  define  FP  to  be  the 
time  integration  of  Q(particle)  over  1  h. 

3. 1 .4.  Q(particle)/Q  T 

[54]  Figure  2d  shows  the  Q(particle)  percentages  of  total 
energy  flux  QT  derived  by  adding  Q(particle)  and  Qy. 
Q(particle)  is  often  much  smaller  than  Qy.  For  cases  of 
significant  Qr  over  timescales  of  several  minutes  or  more, 
the  amount  due  to  Q(particle)  is  generally  less  than  half.  For 
example  on  20  November  there  is  a  single  large  spike  in 
Q(particle)  that  reaches  80  mW/m2.  However,  during  the 
rest  of  the  night  even  though  Q(particle)  is  above  5  mW/m2, 
Qy  is  larger  than  Q(particle)  most  of  the  time. 

3.1.5.  Riometer  Absorption  and  Relative  Energy 

[55]  One  aspect  of  this  study  is  to  determine  whether  the 
storm  on  20  November  shows  a  different  E  than  on  other 
nights  of  this  period.  The  riometer  absorption  (Ab)  data 
shown  in  Figure  2e  does  show  a  significant  difference.  Since 


riometer  absorption  reflects  the  presence  of  high-energy 
particles,  it  is  clear  that  observations  on  20  November  do 
not  resemble  those  on  16  November  or  17  November,  and  in 
fact  20  November  shows  less  absorption  than  on  either  of 
those  days.  The  derivation  of  the  average  electron  energy  is 
shown  in  Figure  2f  using  the  accepted  scaling  of  Ab 
described  above  (E  is  proportional  to  Ab/Q(particle)° 5). 
Since  the  derived  energies  are  only  relative,  the  values  are 
scaled  from  0  to  20  keV.  To  convert  to  absolute  E  values,  a 
scaling  factor  will  be  used  in  later  figures.  On  16  November 
and  17  November  there  are  significant  periods  of  relatively 
high  E  particle  deposition.  This  also  appears  true  on 
18  November  and  19  November.  However,  while  there  are 
certainly  such  periods  on  20  November,  E  appears  consider¬ 
ably  smaller  during  much  of  20  November  and  2 1  November. 

3.1.6.  NOAA  POES  Data 

[56]  Data  from  the  NOAA- 16  and  NOAA- 17  satellites 
allow  a  quick  view  of  the  energetics  of  the  precipitating 
particles.  While  DMSP  satellites  also  provide  these  data,  the 
NOAA- 16  satellite  is  able  to  provide  data  in  a  unique  time 
frame,  namely,  the  postmidnight  sector  near  Poker  Flat. 

[57]  Figure  3  shows  all  of  the  2-s  observations  made  by 
the  NOAA- 16  satellite  in  the  northern  hemisphere  during 
the  period  from  0000  UT  on  16  November  2003  to  2359  UT 
on  21  November  2003  when  the  electron  flux  Q(electron) 
exceeded  1  mW/m2.  Each  observation  is  tagged  by  an  event 
number  that  ranges  from  1  at  the  beginning  of  the  period 
through  just  over  5500  at  the  end  of  the  period. 

[58]  To  provide  an  indication  of  when  these  events 
occurred,  blue  vertical  dotted  lines  mark  the  end  of  each 
day  of  observations  and  the  day  number  is  indicated  just  to 
the  left  of  each  line.  In  addition,  all  odd  orbits  are  shown  in 
black  and  even  orbits  are  shown  in  red,  making  the  events 
associated  with  each  orbit  more  clearly  visible.  Vertical 
dashed  lines  indicate  0800  and  1000  UT  on  20  November 
and  0200  and  0600  UT  on  21  November.  These  are 
interesting  periods  that  will  be  discussed  later  in  the  paper. 

[59]  With  this  approach  it  is  easy  to  visualize  when 
geomagnetic  activity  occurred  because  the  most  active  days 
have  the  most  events.  Note  there  are  many  more  events  on 
20  November  than  on  19  November.  Furthermore,  the 
activity  on  each  orbit  is  indicated  as  the  number  of  events 
per  orbit. 

[60]  Figure  3a  shows  Q(electron)  for  all  the  events.  The 
most  geomagnetically  active  days  16,  17,  and  20  November 
have,  in  general,  higher  fluxes  than  on  the  less  active  days 
1 8  November  and  1 1  November.  The  largest  fluxes  occur  on 
the  superstorm  day  of  20  November. 

[61]  Since  one  focus  of  this  paper  is  on  the  energy  of  the 
precipitating  electrons,  Figure  3b  uses  TED  data  to  show  the 
fraction  of  Q(electron)  due  to  particles  with  energies  below 
1  keV.  While  high  ratios  often  exceeding  0.5  occur  on  all 
nights,  there  is  a  significant  difference  between  the  results 
on  20  November  and  on  the  other  nights.  On  all  other  nights 
than  20  November  the  high-ratio  events  consist  of  one  or 
two  consecutive  observations.  This  suggests  the  presence  of 
narrow  features,  perhaps  auroral  arcs.  However,  after 
~1000  UT  on  20  November  there  are  several  periods  during 
which  this  ratio  exceeded  0.5  over  many  consecutive 
measurements.  This  suggests  a  diffuse-type  enhanced, 
low-energy  aurora  consistent  with  the  presence  of  Type  A 
red  aurora. 
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Figure  3.  Plots  of  the  number  of  2  s  observations  when 
the  electron  flux  is  above  1  mW/m2  for  various  quantities 
from  NOAA  16  data  from  the  period  of  0000  UT  on 
16  November  2003  to  2400  UT  on  21  November  2003 
during  night  northern  hemisphere  overpasses.  Each 
observation  is  tagged  by  an  event  number  that  ranges  from 
1  at  the  beginning  of  the  period  through  just  over  5500  at 
the  end  of  the  period.  To  provide  an  indication  of  when 
these  events  occurred,  blue  vertical  dashed  lines  mark  the 
end  of  each  day  of  observations  and  the  day  number  is 
indicated  just  to  the  left  of  each  line.  In  addition,  all  odd 
orbits  are  shown  in  black  and  even  orbits  are  shown  in  red. 
Thus  the  events  associated  with  each  orbit  are  clearly 
visible.  Vertical  dashed  lines  are  also  shown,  indicating 
0800  and  1000  UT  on  20  November  and  0200  and  0600  UT 
on  21  November.  Shown  are  (a)  Q(electron),  (b)  fraction  of 
O(electron)  below  1  keV,  and  (c)  ratio  of  Q(protons)  to 
Q(total).  Q(total)  is  the  sum  of  Q(electron)  and  Q(proton). 


[62]  Even  more  striking  is  that  on  all  the  previous  nights 
there  are  many  periods  when  the  ratio  goes  nearly  to  0, 
indicating  the  presence  of  few  low-energy  electrons.  The 
period  from  0800  to  1000  UT  on  20  November  shows 
evidence  of  some  enhanced  low-energy  electrons  as  the 
ratio  always  exceeds  zero.  However,  after  1000  UT  the  low- 
energy  component  increases,  resulting  in  a  higher  ratio  well 
above  zero.  In  fact,  from  1000  UT  until  2400  UT  on 
20  November  the  ratio  goes  close  to  zero  only  once.  Thus 
a  persistent  low-energy  component  exists  throughout  this 
period. 

[63]  Figure  3c  shows  the  other  striking  difference  in  the 
auroral  activity  on  20  November.  This  figure  plots  the  ratio 
of  the  proton  flux  to  the  sum  of  the  proton  and  electron  flux. 
The  proton  flux  is  taken  from  the  TED  and  MEPED 
instruments  while  the  electron  flux  is  taken  from  the  TED 
instrument  alone.  Clearly,  there  is  an  enhanced  proton 
component  on  20  November  compared  to  the  other  nights. 
However,  a  careful  examination  of  the  data  does  not  reveal 


any  specific  correlation  between  the  enhanced  proton 
component  and  the  enhanced  low-energy  component.  A 
change  in  the  proton  ratio  is  not  necessarily  correlated  with 
a  low-energy  ratio. 

[64]  Data  from  NOAA- 17,  which  samples  the  northern 
hemisphere  in  the  premidnight  sector  shows  similar 
features.  However,  one  difference  is  that  for  NOAA- 16 
data,  the  period  of  20  November  is  generally  characterized 
by  larger  Q(electron)  values  than  during  the  rest  of  the 
period.  For  the  NOAA- 17  orbit  the  opposite  is  true  and  the 
Q(electron)  values  are  smaller  on  20  November. 

[65]  The  POES  data  taken  together  with  the  riometer 
results  suggest  that  the  large  storm  of  20  November  was 
associated  with  lower-energy  precipitating  electrons.  Alone 
these  data  are  not  able  to  answer  the  questions  associated 
with  composition  change  and  they  do  not  allow  a  continu¬ 
ous  monitor  of  the  electron  energy  distribution.  To  pursue 
this  analysis  requires  investigation  of  the  individual  nights. 

3.2.  Geomagnetic  Activity  Prior  to  the  Storm 

[66]  Data  from  2  d,  17  November  and  19  November,  are 
next  examined  in  detail.  The  first  day  represents  a  disturbed 
period  when  the  DMSP  midnight  auroral  boundary  index 
(MABI)  [Rich  and  Denig ,  1992]  showed  the  southern 
most  portion  of  the  auroral  boundary  to  be  at  a  magnetic 
latitude  of  57°.  On  19  November  the  MABI  indicated  quiet 
conditions  when  the  boundary  was  at  a  magnetic  latitude  of 
60°  to  61°. 

3.2.1.  DMSP  J  Sensor  F16  Data 

[67]  Figure  4  shows  a  series  of  plots  from  the  DMSP  FI  6 
sensor  from  both  17  November  (Figures  4a-4g)  and 
19  November  (Figures  4aa-4gg).  Data  are  restricted  to 
events  when  the  electron  flux  Q(electron)  is  greater  than 
2  mW/m2.  Note  that  unlike  the  POES  data,  here  both 
hemispheres  are  shown.  Figures  4c-4d  and  4cc-4dd  show 
the  latitude  and  longitude  of  each  pass  and  allow  the  reader 
to  identify  northern  and  southern  polar  passes. 

[68]  Because  data  are  presented  as  a  function  of  time, 
individual  points  (plus  signs)  are  closely  spaced  data 
samples  that  to  the  eye  appear  as  nearly  vertical  lines  as 
the  satellite  quickly  passes  through  the  region  of  auroral 
precipitation  in  each  hemisphere.  A  closer  inspection  of 
each  pass  (as  shown  by  the  vertical  distribution  of  points) 
indicates  there  are  usually  two  closely  spaced  lines  repre¬ 
senting  the  two  passages  of  the  auroral  oval  in  each 
hemisphere.  The  average  values  are  shown  over  each 
hemispheric  passage  as  a  large  triangle.  Solid  lines  connect 
the  triangles  to  help  guide  the  reader. 

[69]  On  17  November,  the  nominal  E  (Figure  4a)  is  near 
5  keV  while  most  of  the  flux  values  (Figure  4b)  are  below 
10  mW/m2.  Although  there  is  considerable  spread  in  E,  few 
values  are  at  or  below  1  keV  and  few  values  are  at  or  above 
10  keV.  The  bottom  three  figures  (Figures  4e-4g)  show  the 
fraction  of  the  flux  below  0.1  keV,  between  0.1  and  1  keV, 
and  above  1  keV,  respectively.  For  most  of  these  events,  on 
average,  80  to  95%  of  the  flux  is  above  1  keV,  and  5  to  20% 
of  the  flux  is  between  0.1  and  1  keV.  In  Figure  4e  the  solid 
line  cannot  be  seen  as  the  flux  with  energies  below  0. 1  keV 
is  usually  less  than  1%  of  the  total.  This  is  what  would  be 
expected  for  Gaussian  and  Maxwellian  particle  distributions 
with  low-  and  high-energy  tails  based  on  Strickland  et  al. 
[1993].  For  this  panel  we  plot  a  dashed  line  with  large 
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Figure  4.  Plots  of  various  quantities  from  DMSP  FI 6  from  the  period  of  0400  to  1600  UT  on  (left) 
17  November  2003  and  (right)  19  November  2003  during  all  times  where  Q(electron)  is  above  2  mW/m2. 
Shown  are  (a)  E ,  the  average  energy  of  the  electrons  (b)  Q(electron),  (c)  latitude,  (d)  longitude,  (e)  fraction 
of  flux  below  0. 1  keV,  (f)  the  fraction  of  the  total  flux  between  0. 1  and  1  keV,  and  (g)  fraction  of  total  flux 
above  1  keV.  (aa  to  gg)  Same  as  Figures  4a-4g.  For  all  the  figures  except  Figures  4c,  4d,  4cc,  and  4dd  the 
average  of  the  data  over  each  hemispheric  passage  are  shown  as  triangles  (see  text).  The  solid  line  is 
shown  to  guide  the  reader.  For  Figures  4e  and  4ee  the  diamonds  (connected  by  a  dashed  line)  represent 
that  fraction  of  the  time  during  a  hemisphere  pass  that  the  flux  below  0.1  keV  contributed  more  than  3% 
to  the  total.  Note  that  the  scales  for  the  panels  for  the  same  quantity  are  the  same  from  day  to  day  to 
visually  enhance  any  differences.  See  also  Figure  7. 


diamonds  which  represents  the  fraction  of  the  time  the  flux 
less  than  0.1  keV  is  0.03  (3%)  or  greater  of  the  total.  (Note 
this  plot  is  only  valid  in  the  interior  of  the  diamonds  which 
intersects  the  line.  The  line  is  only  to  guide  the  reader.)  This 
3%  value  was  chosen  because  it  represents  a  typical  mean 
value  found  during  the  20  November  superstorm.  For  the 
night  of  17  November  where  Kp  is  high  this  occurs  only 
10%  or  less  of  the  time. 


[70]  The  data  from  19  November  when  was  low  are 
similar,  but  there  are  some  differences.  Especially  in  the 
earlier  orbits,  there  is  a  tendency  for  the  northern  and 
southern  hemisphere  passes  to  show  different  E  values. 
They  are  much  smaller  during  northern  hemisphere  passes. 
During  these  lower  E  passes  the  flux  greater  than  1  keV 
drops  to  60%  while  the  flux  between  0.1  and  1  keV 
increases  to  40%.  This  is  expected  as  E  during  those  passes 
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Figure  5.  Poker  Flat  photometer  results  from  0600  to  1600  UT  on  17  November  2003.  Derived  results 
are  plotted  when  Q(particle)  is  above  1  mW/m2.  “Quiet”  atmosphere  results  from  Figures  la  and  Id  are 
also  used.  Shown  are  (a)  Q(particle),  (b)  red/blue,  (c)  E  derived  from  red/blue  versus  844.6nm/blue 
photometric  measurements,  (d)  green  solid  line  is  E  derived  from  OI  (557.7  nm)  temperature  (Figure  Id). 
Red  dashed  line  is  E  derived  from  Riometer  absorption  (Ab)  divided  by  Q° 5.  Also  shown  is  (e)  0/N2 
derived  from  red/blue  versus  01(844.6  nm)/blue  (Figure  la.)  The  blue  stars  are  shown  every  3  h  (0600, 
0900,  1200,  1500  UT)  and  are  the  NRLMSIS  0/N2  values.  The  first  star  at  0600  UT  is  at  0.75  and  is 
somewhat  obscured  by  the  y  axis,  (f)  Solid  line  is  F j,  calculated  from  Qj  integrated  over  1  hour. 
Dashed  line  is  ¥P  derived  from  Q(particle)  integrated  over  1  h.  Note  that  the  scales  for  Figures  5a,  5c,  5d, 
and  5e  are  the  same  from  day  to  day  to  visually  enhance  any  differences.  See  Figures  6,  8,  and  9. 


is  below  2  keV.  However,  the  flux  less  than  0.1  keV  still  is 
below  1%  of  the  total  and  the  fraction  of  the  time  this  value 
exceeds  3%  is  even  lower  than  on  17  November. 

3.2.2.  Ground-Based  Poker  Flat  Data 

[7 1  ]  Figures  5  and  6  show  ground-based  data  from  the 
Poker  Flat  observatory  on  these  2  d,  17  November  and  19 
November.  Each  figure  shows  data  for  the  2  d  in  the  same 
format. 

3.2.2.I.  17  November  2003 

[72]  Figure  5a  shows  Q(particle)  derived  from  the  blue 
photometer  data.  The  MSP  data  indicate  weak  proton 
precipitation  throughout  this  night.  The  particle  energy 
flux  is  intense  during  several  short  periods,  reaching  above 
40  mW/m2. 

[73]  The  red/blue  ratio  (Figure  5b)  are  mostly  well  below 
1,  suggesting  moderate  to  high  values  for  E.  This  is 


confirmed  in  Figure  5c  where  the  photometer-derived 
values  (using  red/blue  and  OI  (844.6  nm)/blue  and  the 
“quiet”  atmosphere)  for  E  are  almost  always  at  5  keV  or 
higher,  qualitatively  consistent  with  the  DMSP  data. 

[74]  Figure  5d  shows  two  different  methods  for  deriving 
£,  the  first  uses  the  OI  (557.7  nm)  temperature  method  and 
the  second  uses  the  combination  of  the  riometer 
and  photometer  blue  data.  The  riometer-derived  E  values 
(Ab/Q0  5),  when  arbitrarily  scaled  by  a  factor  of  20,  are  seen 
to  agree  very  well  with  the  temperature-derived  E  values. 
These  latter  results  are  similar  but  slightly  higher  than  the 
results  from  the  photometer  throughout  most  of  the  night. 
Although  all  these  results  show  that  E  is  high  (on  the  order 
5  keV  or  above)  there  is  a  significant  difference  between 
the  results  in  Figure  5c  and  5d  during  the  period  0900  to 
1100  UT.  During  this  time,  a  large  increase  in  E  is  absent 
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Figure  6.  Same  as  Figure  5  except  where  noted  but  for  19  November  2003.  In  Figure  6d  the  dotted  line 
shows  E  calculated  assuming  the  model  atmosphere  T  is  reduced  by  30K.  In  Figure  6e  the  triangles 
(0900,  1200,  and  1500  UT)  represent  the  NRLMSIS  data  scaled  by  0.6.  These  fall  almost  on  top  of,  and 
are  obscured  by,  both  the  data  and  by  the  axis  at  0900  and  1500  UT.  The  stars,  which  represent  the 
unsealed  NRLMSIS  results  are  also  partially  obscured  at  0900  and  1500  UT  by  the  y  axis. 


from  the  data  calculated  using  the  photometric  red/blue 
versus  OI  (8444.6  nm)/blue  technique.  The  reason  for  this 
is  most  probably  the  existence  of  a  significant  high-energy 
tail  component  in  the  energy  distribution  that  is  not  included 
in  the  model.  Recall  that  both  the  riometer  and  temperature 
techniques  are  sensitive  to  higher  energy  particles  while  the 
photometric  technique  is  less  sensitive.  The  reverse  is  true 
for  low-energy  tails.  This  result  shows  the  advantage  of 
having  multiple  measurements  of  the  same  quantity.  It  is 
also  noted  that  for  most  of  the  0600  to  0800  UT  period,  the 
times  where  interference  from  Cassiopeia-A  might  be 
expected,  there  only  appears  to  be  a  small  difference  in 
the  E  values.  However,  the  expected  interference  may  be  the 
cause  of  excess  absorption  that  peaks  between  0630  and 
0700  UT  resulting  in  a  larger  riometer-derived  E  value. 

[75]  Figure  5e  shows  the  0/N2  column  density  derived 
using  the  quiet  atmosphere  plots  of  red/blue  versus 
01(844.6  nm)/blue  from  Figure  la.  This  technique 
should  produce  0/N2  values  reflective  of  the  combined  E 
region  and  F  region  response.  Also  shown  in  Figure  5e  are 
the  0/N2  values  from  NRLMSIS  calculated  every  3  h.  Since 
the  model  uses  3  h  values  in  its  parameterization,  one 


would  not  expect  finer  timescale  variations  to  be  present  in 
the  model  predictions.  However,  since  the  actual  time 
variations  of  0/N2  are  not  well  known,  it  is  of  interest  to 
see  if  the  measured  variations  have  a  shorter  temporal 
period  than  3  h.  The  model  values  do  not  seem  to  agree 
well  in  a  relative  sense  to  the  0/N2  derived  from  the 
photometer  data.  The  model  values  miss  the  decrease  in 
0/N2  seen  near  1200  UT.  Note  though  that  the  data  shown 
in  Figure  5e  are  quite  noisy  with  point  to  point  fluctuations 
often  exceeding  the  less  than  5%  expected  from  statistical 
fluctuations.  In  fact,  examining  the  Q(particle)  data  from 
Figure  5a  shows  many  periods  during  which  the  aurora  is 
quite  dynamic.  As  discussed  earlier  such  periods  can 
produce  spurious  inferences  because  the  photometer  takes 
several  seconds  to  record  all  four  channels  and  auroral 
brightness  changes  during  that  period  will  cause  incorrect 
ratios  to  be  recorded. 

[75]  Finally,  Figure  5f  shows  fluences  Fy  and  FP  from  Q j 
and  Q  (particle),  respectively,  derived  using  the  previous  1  h 
for  an  integration  period.  Note  the  Fj  is  larger  throughout 
most  of  the  period.  The  maximum  occurs  when  0/N2 
minimizes,  consistent  with  previous  studies  [Christensen 
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et  al .,  1997;  Hecht  et  al ,  2006].  For  this  night  the 
total  electromagnetic  energy  fluence,  derived  by  integrating 
Qj  from  0000  to  1600  UT,  is  about  600  J/m2. 

3.2. 2.2.  19  November  2003 

[77]  The  results  for  19  November  in  Figure  6  simply 
reflect  the  lower-energy  deposition  on  this  night.  The  MSP 
data  (not  shown)  again  indicate  weak  proton  precipitation. 
The  red/blue  data  are  generally  below  1  (Figure  6b) 
suggesting  that  the  photometer-derived  E  values 
(Figure  6c)  are  well  above  2  keV,  consistent  with  the  DMSP 
data.  In  Figure  6d  the  riometer-derived  results,  which  use 
the  same  scaling  as  in  Figure  5,  compare  favorably  with  the 
photometer-derived  results  until  about  1300  UT.  However, 
after  1300  UT  the  riometer  results  show  a  significant 
decrease.  The  temperature-derived  results  for  E  (Figure  6d) 
are  higher  than  both  the  riometer  and  photometer  derived 
results.  Reasons  for  this  may  be  that  the  exospheric 
temperature  is  different  than  is  assumed  in  the  model  or 
wave  activity  near  100  km  exists  that  distorts  the  temper¬ 
ature  profile.  To  determine  the  plausibility  of  these 
scenarios,  the  temperature  profile  used  in  Figure  Id  was 
arbitrarily  lowered  by  30  degrees.  The  resultant  E  values 
derived  from  the  OI  (557.7  nm)  temperature  are  shown  as  a 
black  dotted  line  in  Figure  6d.  The  result  is  a  considerable 
reduction  in  the  derived  E ,  similar  to  the  photometric  results 
and  lends  credence  to  the  suspicion  that  the  model  atmo¬ 
sphere  may  not  be  providing  the  correct  temperature 
profile.  Overall,  the  agreement  between  the  temperature 
and  photometer-derived  results  even  after  1300  UT,  and 
their  disagreement  with  the  riometer  results  suggest  an 
absence  of  the  very  high  energy  electrons  responsible  for 
the  riometer  absorption. 

[78]  The  0/N2  column  densities  shown  in  Figure  6e 
indicate  a  small  nightly  variation  with  less  than  a  20% 
difference  between  the  minimum  and  maximum.  As  shown 
by  the  stars,  the  absolute  NRLMSIS  results  are  50%  higher 
and  when  scaled  (shown  as  triangles)  the  relative  variation 
is  quite  similar.  Note  both  that  the  variations  are  on 
the  order  of  10%  suggested  from  statistics  and  that  the 
Q(particle)  data  in  Figure  6a  show  much  less  variation  than 
on  17  November.  The  minimum  in  0/N2  shown  in  Figure  6e 
occurs  just  after  the  peak  of  Fj  shown  in  Figure  6f.  For  this 
night  the  electromagnetic  energy  fluence  Qj  integrated  from 
0000  to  1600  UT  was  about  160  J/m2.  On  the  basis  of 
previous  results  reported  by  Hecht  et  al  [2006]  this  low 
fluence  is  consistent  with  the  small  variation  in  0/N2. 

[79]  The  main  puzzle  is  the  average  absolute  value  of  0/N2, 
which  is  around  0.5  compared  to  the  model  value  of  0.8. 
The  reason  for  this  difference  has  not  been  determined  but 
could  be  related  to  the  extensive  storm-related  composition 
changes  that  occurred  on  previous  nights. 

3.3.  Geomagnetic  Activity  During  the  Storm 

[80]  Next  we  present  similar  results  for  the  two  storm 
days,  20  November  and  21  November.  On  20  November  the 
MABI  indicates  a  major  storm  with  a  boundary  ranging 
from  a  magnetic  latitude  of  65°  prior  to  the  storm  onset 
around  0800  UT  and  dropping  to  about  48°  at  the  peak  of 
the  storm.  On  21  November  the  MABI  indicates  disturbed 
to  quiet  conditions  with  the  boundary  retreating  from  a 
magnetic  latitude  of  55°  at  0000  UT  to  61°  later  in  the  day. 


3.3.1.  DMSP  J  Sensor  F16  Data 

[si]  Figure  7  shows  J  sensor  data  for  these  2  d.  The  data 
are  plotted  in  a  similar  manner  as  in  Figure  4.  We  remind 
the  reader  that  the  triangles  in  Figures  4e  and  4ee  indicate 
the  fraction  of  the  time  the  flux  below  0. 1  keV  is  0.03  or 
greater  of  the  total  Q. 

[82]  On  20  November  there  is  considerable  flux 
(Figure  4b)  often  exceeding  10  mW/m2  throughout  the 
day  both  before  and  after  the  storm  onset.  In  fact  for  the 
pass  through  the  auroral  zones  near  0813  UT,  the  flux 
reaches  75  mW/m2  and  may  represent  the  shock  discussed 
in  the  next  section.  The  value  of  the  average  electron 
energy,  E  (Figure  4a),  during  the  first  four  passes  through 
the  auroral  zone  are  quite  different  between  the  northern  and 
southern  hemisphere.  This  is  probably  because  the  satellite 
samples  local  time  sectors  that  are  different  in  the  north 
from  the  south.  For  the  more  energetic  passes  the  averages 
are  near  5  keV,  similar  to  what  was  observed  on 
17  November.  For  the  first  and  third  more  energetic  passes 
before  0800  UT  the  fraction  of  the  flux  in  each  energy  band 
(Figures  4e,  4f,  and  4g)  looks  similar  to  the  previous  days, 
with  little  flux  below  1  keV  (less  than  0.01).  The  fifth  pass 
near  0800  UT  is  also  similar,  although  the  flux  below  1  keV 
has  slightly  increased. 

[83]  Beginning  at  1000  UT  the  E  values  drop  with  large 
fractions  occurring  below  1  keV.  Except  for  the  periods 
around  1400  and  1600  UT,  the  average  E  is  at  or  below 
2  keV  with  few  if  any  points  above  5  keV.  However,  these 
changes  are  not  just  because  the  fraction  of  the  flux  between 
0. 1  and  1  keV  has  increased  while  the  fraction  above  1  keV 
has  decreased  as  was  found  for  19  November.  Now 
fractions  for  the  flux  below  0. 1  keV  are  several  percent  of 
the  total,  sometimes  reaching  50%.  The  fraction  of  the  time 
that  the  flux  below  0.1  keV  is  0.03  or  greater  is  typically 
from  0.2  to  as  high  as  0.8. 

[84]  These  results  show  that  the  energy  distribution  is 
clearly  different  than  on  the  previous  nights  for  most  of  the 
period  from  1000  to  2400  UT.  There  is  not  just  a  lower  E 
but  an  enhancement  in  the  low-energy  portion  of  the 
spectrum  below  0. 1  ke V.  There  is  some  variation,  however, 
in  the  energetics.  The  period  around  1400  and  1600  UT 
shows  more  energetic  particles. 

[85]  The  data  from  21  November  (Figures  4aa  to  4gg) 
show  some  similarity  to  that  on  20  November.  During  the 
period  from  0000  to  0600  UT,  the  average  E  values  are  low, 
around  2  keV  or  less.  However,  the  fraction  of  flux  below 
0.1  keV  (Figure  4ee)  does  not  resemble  the  period  from 
1000  to  2400  UT  on  20  November.  Those  fractions  are 
somewhat  enhanced  compared  to  19  November  but  well 
below  those  seen  on  20  November.  After  about  0800  UT 
they  resemble  those  found  on  19  November. 

3.3.2.  Ground-Based  Poker  Flat  Data 
3.3.2.I.  Results  for  20  November  2003 

[86]  The  results  shown  in  Figure  8  are  derived  using  the 
“disturbed”  atmosphere,  except  where  noted.  Figure  8a 
shows  a  large  impulse  of  auroral  energy  at  ~08 10  UT.  This 
is  presumably  related  to  the  shock  observed  by  SOHO  at  LI 
at  0728  UT  as  reported  by  Gopalswamy  et  al  [2005b].  The 
transit  time  of  the  shock  to  the  top  of  the  magnetopause  is 
between  30  and  60  min.  After  0800  UT  the  particle  flux 
decreases  to  a  minimum  around  1030  UT.  There  are  several 
periods  after  this  of  intense  auroral  activity. 
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Figure  7.  Same  as  Figure  4  except  as  noted.  Plots  of  various  quantities  from  DMSP  F16  from  the  period 
of  (a  to  g)  0400  to  2400  UT  on  20  November  2003  and  (aa  to  gg)  0000  to  1 600  UT  on  2 1  November  2003 
during  all  times  where  Q(electron)  is  above  2  mW/m2.  The  triangles  and  diamonds  are  described  in 
Figure  4.  Note  that  the  scales  for  the  same  quantity  are  the  same  from  day  to  day  to  visually  enhance  any 
differences  especially  with  the  results  shown  in  Figure  4. 


[87]  There  appear  to  be  two  regimes  for  the  energetics  on 
this  night.  Prior  to  1000  UT  (and  including  the  period 
of  large  particle  deposition  around  0800  to  0900  UT)  the 
red/blue  ratio  in  Figure  8b  is  low,  indicative  of  moderate  to 
high  values  of  E.  These  are  consistent  with  those  DMSP 
results  which  sample  the  nighttime  auroral  oval  which  for 
this  period  occurs  in  the  southern  hemisphere.  Just  before 
1000  UT  this  ratio  increases  dramatically  to  values  of  3  to 
4  and  it  remains  at  that  level  or  higher  for  the  remainder  of 
the  night.  The  peak  in  the  red/blue  ratio  is  15,  although  this 
peak  occurs  during  a  period  of  very  low  flux  suggestive  of 


polar  rain.  However,  at  other  times  (e.g.,  1300  UT)  the 
increase  in  this  ratio  occurs  coincident  with  a  significant 
flux  of  particles.  The  absolute  red  intensity  exceeds  10  kR 
during  these  times.  These  are  all  characteristics  of  great  red 
aurorae  [Vallance  Jones,  1992]. 

[88]  Figure  8b  also  shows  the  relative  proton  flux  given 
by  the  ratio  of  H-Beta  to  blue,  which  has  been  scaled 
upward  by  a  factor  of  20  for  clarity.  Certainly  proton 
precipitation  is  enhanced  after  0800  UT  through  the 
rest  of  the  night.  However,  there  is  no  clear  correlation 
with  red/blue.  There  are  periods  when  both  increase,  sugges- 
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Figure  8.  Same  as  Figure  5  except  where  noted  but  for  20  November  2003.  The  “disturbed” 
atmosphere  results  are  used.  In  Figure  8b  the  green  line  represents  the  ratio  of  H-Beta  to  blue  from  MSP 
data  scaled  upward  by  a  factor  of  20.  In  Figures  8c  and  8e  the  red  dotted  line  is  for  the  “quiet” 
atmosphere. 


ti  ve  of  the  presence  of  enhanced  fluxes  of  low-energy  electrons 
and  protons.  There  are  other  periods  where  no  such  correlation 
exists.  Note  especially  that  the  proton  flux  enhancement 
between  0800  and  1 000  UT  does  not  correlate  with  the  changes 
in  red/blue  seen  during  that  same  period. 

[89]  The  values  of  E  derived  from  the  three  techniques 
(Figures  8c- 8d)  show  some  interesting  trends.  E  derived 
from  the  red/blue  versus  OI  (844.6  nm)/blue  photometer 
technique  (Figure  8c)  shows  values  of  5  to  10  keV  prior  to 
0930  UT,  dropping  to  1  to  2  keV  after  1000  UT  and  the  rest 
of  the  night.  This  seems  consistent  with  the  DMSP  results. 
The  temperature-derived  E  values  shown  in  Figure  8d  also 
show  a  decline,  but  not  nearly  as  large  with  few  values 
below  5  keV.  This  is  not  surprising  because  the  temperature 
technique,  as  noted  earlier,  is  relatively  insensitive  to  the 
low-energy  electrons,  which  would  not  strongly  excite  the 
green  emission. 

[90]  Also  plotted  in  Figure  8c  are  the  E  results  derived 
using  the  “quiet”  atmosphere.  There  are  small  differences 
with  the  energy(0/N2)  results  in  the  quiet  atmosphere  being 
between  about  10(5)%  lower  suggesting  that  the  results 


using  the  photometric  technique  are  relatively  independent 
of  model  atmosphere. 

[91]  The  riometer-derived  E  values  are  also  shown  in 
Figure  8d.  Since  the  riometer  data  are  not  sensitive  to  low- 
energy  electrons,  any  such  CNA  should  be  interpreted  as  an 
indicator  of  high-energy  electrons.  Over  a  significant 
portion  of  the  night,  the  riometer-derived  E  values  are 
considerably  lower,  suggesting  a  relative  absence  of 
the  high-energy  electrons  needed  to  produce  CNA.  Inter¬ 
estingly,  in  the  period  from  0600  to  0700  UT,  high-energy 
electrons  appear  and  then,  just  before  the  storm  onset,  this 
high-energy  component  disappears.  Examination  of  all  the 
imaging  riometer  pixels  (not  just  those  in  the  magnetic 
zenith)  strongly  suggest  that  because  the  enhancement  in 
CNA  prior  to  0700  UT  is  an  artifact,  due  to  the  passage  of 
the  radio  star  Cassiopeia-A.  Ignoring  the  artifact,  the  period 
when  high-energy  electrons  are  present  occurs  almost 
impulsively  around  0800  UT  when  the  shock  arrives. 
Some  high-energy  electrons  remain  after  this  event,  but 
after  1000  UT  they  diminish  considerably  until  ~1200  UT. 
The  poor  agreement  between  the  riometer-derived  E  values 


17  of  28 


A01310 


HECHT  ET  AL.:  THERMOSPHERIC  COMPOSITION  ON  17-21  NOV  2003 


A01310 


40 

%  30 

a 

3.5 
c;  2.8 

m 

^  20 

- 

5  2.1 

1  .  .  - 

i  10 
0 

_ : 

TJ  1.4 

V 

“  0.7 
0.0 

4  5  6  7  8  9  10  11  12  13  14  15  16 


> 

c 

j* 

w 


20 

15 

10 

5 

0 


4  5  6  7  8  9  10  11  12  13  14  15  16 
Green  Solid  (T),  Red  Dash  (Rio) 


d 

l 

\  r  Mr 

ink . ■ 

: 

1  Lf  v 

4  5  6  7  8  9  10  11  12  13  14  15  16 


11/21  UT  Hour 


Solid(Joule),  Dash(Particle) 


r*  50 

£  40 

f 

ft  'V 

- 

e(J 

CO 

0 

t 

U  U 

- 

c  20 
0) 

/ 

-  / 

/ 

✓  J 

It 

1  \ 

Flu 

O  O 

mil 

/ 

4  5  6  7  8  9  10  11  12  13  14  15  16 


11/21  UT  Hour 


Figure  9.  Same  as  Figure  8  but  for  0400  to  1600  UT  on  21  November  2003.  The  “disturbed” 
atmosphere  results  are  used. 


and  the  other  techniques  between  0700  and  0800  UT 
suggest  an  absence  of  very  high-energy  electrons  during 
this  period.  The  better  agreement  of  the  riometer- 
derived  and  temperature-derived  E  values  after  1200  UT 
suggest  the  presence  of  a  more  nominal  high-energy 
electron  component. 

[92]  The  0/N2  results,  which  are  derived  from  the 
photometer  data,  are  shown  in  Figure  8e.  Between  0600 
and  0700  UT,  Figure  8e  shows  values  just  above  0.6,  a  value 
slightly  higher  than  was  found  at  the  end  of  19  November. 
This  value  gradually  increases  to  0.75  and  then  decreases 
almost  steadily  until  1400  UT  to  about  0.2.  It  then  increases  to 
about  0.25  at  1 600  UT.  Also  shown  are  the  NRLMSIS  values 
which  start  at  about  0.9  at  0600  UT.  These  values  decline  to 
0.75  at  0900  UT  and  continue  to  decline  reaching  a  minimum 
of  0.45  at  1 5  UT.  Similar  to  what  was  found  for  E,  the  results 
using  the  “quiet”  atmosphere  with  the  enhanced  Gaussian 
are  between  5  and  10  percent  lower  than  when  the 
“disturbed”  atmosphere  is  used. 

[93]  There  are  several  aspects  of  these  results  worth 
noting.  First,  the  decline  of  the  photometric-derived  0/N2 
results  is  generally  smooth.  It  does  not  show  any  jump  at 
1000  UT  associated  with  the  change  in  particle  energetics. 
Second,  NRLMSIS  shows  about  a  factor  of  two  decrease 


while  the  data  show  about  a  factor  of  3  decrease.  This  is 
comparable  to  what  was  found  during  the  great  magnetic 
storm  of  February  1986  [Hecht  et  al ,  1991].  Third,  these 
results  are  reassuring  in  that  the  composition  results  and  the 
E  values  derived  from  the  photometer  technique  (Figures  8c 
and  8e)  are  fairly  insensitive  to  the  choice  of  model 
atmosphere.  Such  a  comparison  was  difficult  to  perform 
previously  when  the  fG  scaling  factor  was  used  [Hecht  et  al ., 
1991]. 

[94]  Finally,  Figure  8f  shows  F j  and  F/>  Again  F j  is 
significantly  higher  through  most  of  the  period.  The  F j 
values,  even  when  low  (after  0800  UT),  are  generally  much 
larger  than  on  previous  nights  and  are  consistent  with  a 
continuous  decline  in  0/N2  seen  throughout  the  night.  The 
smallest  F,  values  (after  0800  UT)  occur  around  1 500  UT 
and  could  be  the  reason  for  the  slight  increase  in  0/N2 
around  that  time.  The  electromagnetic  energy  fluence 
obtained  by  integrating  Qy  over  the  0000  to  1600  UT  period 
is  2700  J/m2,  the  largest  value  seen  of  all  the  observations 
reported  on  here  and  by  Hecht  et  al  [2006]. 

3.3.2.2.  Results  for  21  November  2003 

[95]  Figure  9  shows  the  results  for  21  November  which 
covers  the  period  after  the  peak  of  the  storm.  As  seen  in 
Figure  9a,  only  one  large  storm-related  auroral  event  occurs 
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Figure  10.  Plots  of  Fort  Yukon  photometer  results.  Except 
for  Figure  lOd,  first  point  is  at  0600  UT.  (a)  0/N2  column 
density  ratio  on  17  November  2003.  (b)  Same  as  Figure  10a 
on  19  November  2003.  (c)  Same  as  Figure  10a  on 
20  November  2003.  (d)  E  on  20  November  2003.  (e)  Same 
as  Figure  10a  on  21  November  2003. 


and  that  period  is  between  0500  and  0700  UT.  During  this 
time  the  red/blue  ratio  shown  in  Figure  9b  is  lower  than  on 
20  November  with  values  around  1  or  less.  The  E  values 
shown  in  Figure  9c  are  higher  than  on  20  November. 
Between  0400  and  0500  UT  the  red/blue  ratio  declines 
from  values  above  2  to  less  than  1 .  The  temperature-derived 
E  values  shown  in  Figure  9d  are  in  reasonable  agreement 


with  the  photometer  derived  values  in  Figure  9c.  The 
riometer-derived  E  values  also  shown  in  Figure  9d  indicate 
an  absence  of  high-energy  particles  prior  to  0500  UT. 

[96]  The  0/N2  ratios  both  from  the  NRLMSIS  model  and 
from  the  data  shown  in  Figure  9e  are  only  slightly  above 
those  found  on  20  November  in  the  post- 1000  UT  period. 
The  Fy  values  shown  in  Figure  9f  are  closer  to  that  seen 
on  19  November  than  to  the  other  nights.  The  total 
electromagnetic  energy  fluence  for  the  night  obtained  from 
integrating  Q j  from  0000  to  1600  UT  is  about  120  J/m2,  the 
lowest  of  the  four  nights  presented  here.  Despite  this  small 
amount  of  heating,  the  0/N2  ratios  do  not  recover  to  their 
presuperstorm  levels. 

3.4.  Other  Data 

3.4.1.  Satellite  Particle  Data  on  20-21  November 

3.4.1. 1.  NO  A  A  POES  Data 

[97]  The  POES  data  from  NOAA-16  were  shown  in 
Figure  3  and  discussed  earlier.  They  show  a  marked 
increase  in  low-energy  particles  after  1000  UT  until 
2400  UT  on  20  November.  There  is  also  a  minimum  in 
the  low-energy  component  around  1600  UT  consistent  with 
the  DMSP  FI 6  data.  On  21  November  the  data  show 
increased  low-energy  particles  between  about  0200  and 
0600  UT. 

[98]  The  data  for  NOAA  17  (not  shown)  are  similar  to  the 
NOAA-16  results.  The  only  difference  in  these  data  is  that 
the  enhancement  of  the  low-energy  electrons  after  2000  UT 
is  larger  than  in  the  NOAA-16  results. 

3.4.1.2.  DMSPF15 

[99]  The  data  from  the  DMSP  FI 5  J  sensor  also  show  an 
increase  during  20  November  in  the  very  low  energy 
(<0.1  keV)  electrons  but  not  to  the  same  magnitude  as 
was  seen  in  FI 6.  While  the  increase  begins  at  0800  UT,  the 
fraction  of  the  time  these  low-energy  electrons  are  above 
0.03  and  thus  carry  more  than  3  percent  of  the  total  flux, 
reaches  a  local  maximum  of  around  0.1  at  approximately 
1030  UT.  It  stays  around  that  value  for  the  remainder  of 
20  November.  There  are  some  periods  where  this  value 
(the  fraction  of  time  the  low-energy  particles  carry  more 
than  3  percent  of  the  total  flux)  is  higher,  at  almost  0.4,  and 
lower,  close  to  0.01  near  1400  UT.  However,  in  general  this 
value  is  enhanced  over  the  pre-0800  UT  period,  but  below 
that  seen  in  the  FI 6  data.  On  21  November  (except  for  one 
orbit  near  0200  UT  that  value  is  also  enhanced  in  the  FI 6 
data)  this  value  is  at  levels  similar  to  the  pre-0800  UT 
values  on  20  November. 

[100]  The  J  sensor,  however,  does  not  appear  to  record  the 
large  flux  seen  on  the  ground  in  the  Alaska  data  and  on 
DMSP  FI 6.  This  is  probably  because  the  pass  through  the 
auroral  zone  occurs  near  0804  UT,  which  is  before  the  large 
increase  was  recorded  at  Poker  Flat. 

3.4.2.  Ground-Based  Data 
3.4.2.I.  Fort  Yukon  Photometer 

[101]  Selected  results  from  this  instrument  are  shown  for 
all  four  nights  in  Figure  10.  The  composition  results  shown 
in  Figures  10a- 10c  and  Figure  lOe  are  similar  to,  but 
not  identical  with,  those  seen  at  Poker  Flat.  This  is  not 
surprising  as  the  Joule  heating,  which  appears  to  control 
most  of  this  change,  has  a  slightly  different  time  history  at 
the  two  sites.  However,  that  the  decrease  in  0/N2  is  largest 
on  20  November  is  confirmed  by  these  data.  Of  particular 
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Figure  11.  Plots  of  various  quantities  from  (left)  Gillam  and  (right)  Pinnawa  from  meridian  scanning 
photometer  data  during  the  period  of  0000  to  1200  UT  on  20  November  2003  where  Q  (particle)  is  above 
1  mW/m2,  showing  (a)  Q(particle)  derived  from  N2+  data,  (b)  Q(proton)/(Q(particle),  (c)  red/blue,  and  (d) 
E  derived  using  OI  (557.7  nm)/blue  versus  red/blue  using  Figure  lc. 


interest  are  the  E  data  on  20  November  shown  in  Figure  lOd. 
As  is  found  at  Poker  Flat,  there  is  a  significant  decrease  in  E 
after  about  0930  UT,  and  low  values  between  1  and  3  keV 
are  measured  throughout  the  rest  of  the  night. 

3.4.2.2.  Pinnawa  and  Gillam  Photometers 

[102]  To  further  investigate  the  large-scale  nature  of  this 
low-energy  particle  precipitation,  Figure  11  shows  results 
from  the  two  mid-Canadian  photometer  sites  at  Pinnawa 
and  Gillam.  Figures  11a  and  llaa  show  the  particle  flux 
Q(particle)  derived  from  the  N2  data.  Also  shown  are  the 
proton  fraction  (Figures  lib  and  llbb),  the  red/blue  ratio 
(Figures  11c  and  llcc),  and  E  (Figures  lid  and  lldd). 

[103]  As  is  found  at  Poker  Flat  and  at  Fort  Yukon,  both 
Canadian  sites  show  increases  in  red/blue  and  decreases  in 
E  after  0900  to  1000  UT.  Around  0900  UT  the  E  value  is 
5  keV,  which  decreases  to  close  to  1  keV  after  1000  UT. 
This  is  consistent  with  what  was  found  at  the  two  Alaska 
sites. 

[104]  One  additional  aspect  of  these  data  is  also  worth 
discussing.  Figures  lib  and  llbb  show  the  ratio  of  the 
proton  flux  derived  from  H-beta,  with  the  total  flux  derived 


from  the  N2  photometer.  Both  use  constant  values  to 
convert  the  measured  brightness  to  flux.  In  the  case  of  the 
N2  data  these  values  are  somewhat  insensitive  to  E  or  to 
protons  or  electrons  [Strickland  et  al.y  1993],  as  discussed 
earlier.  However,  the  H-beta  conversion  is  sensitive  to  the 
proton  average  energy  and  this  uncertainty  probably 
accounts  for  the  values  above  1  seen  in  Figure  lib. 

3.4.3.  GUVI  Data 

[105]  As  noted  earlier,  an  extensive  analysis  of  GUVI 
dayglow  data  has  been  reported  for  this  period  [Meier  etal., 
2005;  Crowley  etal. ,  2006].  Figure  1 2  shows  the  day  side  0/N2 
column  density  ratio  for  20  November.  Although  these  orbits 
cover  the  entire  24  h  UT  period,  the  local  solar  time  for  the 
TIMED  satellite  during  these  dayside  periods  are  generally 
between  about  1000  and  1200  LST.  (Specifically,  a  nominal 
southbound  equatorial  crossing  occurs  about  1150  LST.  At 
40N,  50N,  and  60N  these  crossings  occur  around  1105, 1040, 
and  1000  LST,  respectively.)  Prior  to  1000  UT  the  ratio  has  a 
value  of  about  1,  consistent  with  geomagnetically  quiet 
conditions.  However,  just  after  1000  UT  a  large  decrease  in 
0/N2  appears,  propagating  from  the  north  toward  the  equator  at 


20  of  28 


A01310 


HECHT  ET  AL.:  THERMOSPHERIC  COMPOSITION  ON  17-21  NOV  2003 


A01310 


Day  324  (11/20)  2003 
GUV1  O/N: 
x 


•180  -90  0  90  X  180 


UT  19.49  16.26  13.02  9.78  6.55  3.31  0.07 


Figure  12.  Dayside  0/N2  column  ratio  on  20  November 
2003  derived  from  GUVI  data.  The  two  points  marked  as  a 
red  cross  are  the  north  and  south  magnetic  poles.  The  mean 
UT  times  for  equatorial  crossings  are  shown  on  the  black 
scale  below. 


later  UT.  It  is  of  interest  to  note  that  no  such  decrease  is  apparent 
in  the  period  from  0800  to  1 000  UT. 

[106]  Two  orbits  of  nightside  auroral  zone  GUVI  data 
are  most  relevant  to  this  study  and  these  are  shown  in 
Figure  13.  Figures  13a  and  13b  show  plots  of  GUVI  derived 
results  for  E  and  Q  on  20  November  from  orbit  10559, 
which  passes  over  Western  Canada  at  about  0915  UT 
and  from  the  next  orbit  which  passes  over  Alaska  around 
1055  UT.  The  black  dots  in  Figures  13c  and  13d  locate  the 
Poker  Flat  and  Fort  Yukon  ground  sites.  The  white  arrows 
show  the  satellite  direction.  The  Q  plots  show  the  extent  of 
the  aurora  which  just  reaches  the  US-Canadian  border,  a 
somewhat  lower  latitude  than  would  be  expected  from  the 
DMSP  MABI.  The  E  plots  show  a  considerable  reduction  in 
E  between  these  two  orbits.  This  is  qualitatively  consistent 
with  the  other  data  sets.  However,  since  the  GUVI 
algorithm  cannot  reliably  obtain  average  energies  when  they 
are  below  2  keV,  the  magnitude  of  E  more  resembles  that  seen 
from  the  temperature-derived  E  values  in  Figure  8d.  Unlike 
the  single  point  photometer  measurements,  the  GUVI  data 
can  be  used  to  determine  the  geographical  extent  over  which 
the  low-energy  particles  were  present. 

4.  Discussion 

4.1.  Composition  Changes  During  This  Storm  Period 

[107]  Great  magnetic  storms  are  rare  events  and  prior  to 
this  study,  the  only  composition  change  data  observed  was 
during  the  1986  storm  [Hecht  et  al ,  1991].  The  present 
storm  is  the  first  one  since  the  four-channel  photometric 
technique  was  adopted. 
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Figure  13.  GUVI  images  of  derived  data  on  20  November  for  two  orbits,  (left)  10559  at  approximately 
0915  UT  and  (right)  10560  at  approximately  1055  UT.  Shown  is  (a)  E  for  10559,  (b)  Q  for  10559,  (c)  E 
for  10560,  and  (d)  Q  for  10560.  The  black  dots  in  Figures  13c  and  13d  locate  the  Poker  Flat  and  Fort 
Yukon  sites.  The  white  arrows,  seen  in  the  left  comer  of  each  swath,  show  the  satellite  direction. 
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4.1.1  Model  Versus  Data  Comparison 

[  108]  The  NRLMSIS  model  results  show  small  changes  in 
0/N2  during  19  November,  slightly  larger  changes  on 
17  November,  and  the  largest  changes  on  20  November. 
These  trends  follow  the  Kp  and  F j  values  for  these  days.  On 

19  November  the  model  shows  that  the  0/N2  values  had 
recovered  from  the  storm  activity  on  17  November.  The 
model  indicates  a  significant  decline  in  the  0/N2  values  on 

20  November  that  only  slightly  recovered  to  prestorm 
values  on  21  November. 

[109]  The  data  show  a  somewhat  more  dynamic  picture. 
On  17  November,  the  absolute  value  of  0/N2  agreed  with 
the  model  prior  to  the  significant  activity,  but  there  was  a 
much  greater  depletion  in  0/N2  than  seen  by  the  model.  On 
19  November  the  inferred  changes  in  0/N2  are  small  and 
comparable  to  the  model  predictions  although  the  absolute 
values  are  below  the  model.  This  suggests  that  recovery  had 
not  completely  occurred  from  the  prior  day’s  geomagnetic 
activity.  On  20  November  both  model  and  data  show 
continuous  decline  during  the  night,  although  the  data  show 
a  larger  decrease  of  about  a  factor  of  three  compared  to 
the  factor  of  two  predicted  by  the  model.  The  data  on 

21  November  show  that  0/N2  had  not  recovered  from  the 
superstorm  activity  on  20  November,  a  result  replicated  by 
the  model. 

[no]  Interestingly,  the  model  seems  to  do  a  reasonable  job  of 
reproducing  the  changes  in  0/N2  when  there  is  little  geomag¬ 
netic  activity  and  during  larger,  more  sustained  events  such  as 
on  20  November.  However,  the  model  does  a  poorer  job  during 
moderate  storms  such  as  on  1 7  November.  It  appears  as  if  the 
model,  which  uses  3  h  a p  values,  does  not  capture  large 
geomagnetic  effects  on  short  timescales.  Even  during  large 
storms  such  as  on  20  November  there  is  the  tendency  to 
underpredict  the  changes  in  0/N2. 

4.1.2.  Causes  of  Composition  Change 

[m]  As  Christensen  et  al.  [1997]  have  pointed  out,  the 
qualitative  effects  of  auroral  energy  input  have  been  known 
for  a  long  time.  Energy  is  deposited  into  the  atmosphere  due 
to  dissipation  of  electric  currents,  known  as  Joule  heating, 
and  due  to  the  dissipation  of  kinetic  energy  from  precipi¬ 
tating  particles.  These  energy  losses  cause  circulation  cells 
to  form,  vertical  winds  to  increase  in  the  vicinity  of  the 
heating,  and  subsequent  uplifting  of  air  parcels  containing 
relatively  little  atomic  oxygen.  At  a  given  altitude  above  the 
heating  this  results  in  a  reduction  of  the  ratio  of  light 
molecular  weight  species  to  heavy  molecular  weight  species 
[e.g.,  Hays  et  al .,  1973;  Prolss ,  1980].  Certainly,  for  large 
geomagnetic  storms  and  at  altitudes  above  150  km,  such 
effects  are  found  in  the  predictions  of  the  MSIS  models. 
Similar  effects  should  also  be  seen  on  a  more  local  scale 
during  and  after  auroral  substorms.  Thus  deposition  of 
auroral  energy  should  drive  vertical  circulation,  resulting 
in  changes  in  atmospheric  composition.  What  is  not  known 
is  how  well  current  models  reproduce  actual  composition 
change  during  auroral  activity. 

[112]  On  the  3  d  with  significant  declines  in  0/N2  during 
the  night,  F j  greatly  exceeds  F^  Thus  it  appears  that  one 
major  cause  of  changes  in  0/N2  is  the  magnitude  of  F j. 
Clearly  transport  effects  may  also  be  important  [Hecht  et  al , 
2000]  and  that  may  be  partly  responsible  for  the  low  0/N2 
values  on  21  November  and  on  19  November,  but  more 


work  needs  to  be  done  to  quantify  the  importance  of 
transport. 

4.2.  Magnetic  Cloud  and  Low-Energy  Particles 

[113]  Sivjee  and  Shen  [1997]  proposed  that  the  precipita¬ 
tion  of  low-energy  particles  whose  average  energies  are 
considerably  below  1  keV  may  be  associated  with  the 
passage  of  magnetic  clouds.  Those  results  were  based  on 
temperatures  inferred  from  optical  emissions  observed  over 
Kangerlussuaq,  Greenland  (67N,  209E  geographic,  74.5N, 
41.6W  geomagnetic).  Like  the  present  observations,  that 
event  occurred  during  a  major  magnetic  storm  in  which  Dst 
reached  —  1 30  nT.  The  new  data  presented  for  the  storm  on 
20  November  2003  support  the  thesis  of  Sivjee  and  Shen 
[1997]  that  there  is  an  enhancement  of  low-energy  particles 
during  a  magnetic  cloud  event. 

[114]  The  plasma  properties  in  the  solar  wind  can  be 
determined  from  spacecraft  data.  Gopalswamy  et  al. 
[2005b]  notes  that  the  initial  solar  wind  shock  was  present 
at  LI  around  0728  UT.  Allowing  approximately  0.7  h  to 
advance  the  shock  from  LI  to  the  magnetopause,  the  period 
just  after  0810  UT  should  correspond  to  observations  of 
enhanced  auroral  electron  flux.  This  is  indeed  what  was 
seen  by  the  ground-based  photometers  at  Poker  Flat  at 
0810  UT  and  the  DMSP  F16  J  sensor  around  0813  UT. 
However,  the  F15  J  sensor  that  passed  through  the  auroral 
zone  around  0804  UT  did  not  see  an  enhanced  auroral 
electron  flux  since  that  probably  occurred  before  the  shock 
reached  the  magnetopause. 

[115]  Gopalswamy  et  al.  [2005b]  discusses  a  magnetic 
cloud  associated  with  an  interplanetary  coronal  mass  ejec¬ 
tion  (ICME)  that  was  present  around  the  Earth  starting  at 
approximately  1000  to  1100  UT  on  20  November  2003  and 
ending  around  0000  UT  on  21  November  2003.  The 
passage  of  the  cloud  was  marked  by  enhanced,  smoothly 
rotating  magnetic  fields,  enhanced  proton  density,  and  a 
greatly  depressed  plasma  /?,  which  is  the  ratio  of  the  bulk 
plasma  thermal  pressure  to  the  magnetic  field  pressure, 
relative  to  the  ambient  solar  wind  [Gopalswamy  et  al , 
2005b,  Figure  1].  These  three  features  characterize  magnetic 
clouds  [Burlaga,  1991]. 

[1 16]  Although  the  magnetic  cloud  comprised  most  of  the 
disturbed  interval  inside  the  ICME,  there  were  two  addi¬ 
tional  periods  of  enhanced  proton  density  and  reduced 
plasma  (3.  The  first  occurred  at  the  end  of  the  sheath  region 
between  the  shock  and  the  cloud  where,  between  0900  and 
1000  UT,  the  plasma  (3  is  sharply  reduced.  The  second  is 
following  the  cloud  and  extends  from  approximately 
0000  UT  to  0600  UT  on  21  November  2003.  Although 
the  plasma  (3  remains  below  1  throughout  this  interval,  it 
should  be  noted  that  this  period  is  associated  with  an 
enhancement  in  the  plasma  (3  compared  to  the  magnetic 
cloud  interval.  H-alpha  images  from  SOHO  during  the  solar 
event  show  a  dense  filamentary  structure  ejected  along  with 
the  CME.  The  subsequent  observation  of  the  enhanced 
density  region  in  the  ACE  data  is  consistent  with  such 
filamentary  material,  which  is  often  found  in  the  trailing 
edge  of  ICMEs  [ Gopalswamy  et  al .,  1998]. 

[117]  As  noted  above  most  of  the  data  presented  in  this 
paper  support  the  assertion  by  Sivjee  and  Shen  [1997]  that 
the  precipitation  of  lower  average  energy  particles  may  be 
associated  with  magnetic  clouds.  The  ground-based  data 
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from  both  Alaska  and  Canadian  sites  show  a  decrease  in  E 
in  the  period  from  0900  to  1000  UT.  This  decrease  in  E 
persists  through  the  rest  of  20  November.  The  Alaska  data 
show  low  E  values  until  about  0400  UT  on  21  November, 
after  which  they  begin  to  increase.  Taken  alone,  these  data 
would  suggest  that  lower  E  values  are  associated  in  some 
way  with  the  passage  of  the  core  of  the  magnetic  cloud 
occurring  from  1000  to  2400  UT  on  20  November.  How¬ 
ever,  looking  at  the  boundary  regions  of  the  cloud,  the  data 
suggest  a  much  sharper  boundary  for  the  low-energy 
particles  associated  with  the  front  of  the  cloud  and  a  more 
diffuse  boundary  on  the  backside  of  the  cloud.  The  NOAA 
data  from  Figure  3b,  which  show  a  persistent  enhanced  low- 
energy  electron  flux  component  during  the  passage  of  the 
magnetic  cloud,  also  support  this  picture. 

[i is]  The  DMSP  J  data  and  the  red/blue  data,  however, 
suggest  that  this  picture  is  incomplete.  The  enhancement  of 
the  red/blue,  although  an  indicator  of  a  lower  E  value,  may 
also  be  indicative  of  enhanced  precipitation  of  soft  electrons 
with  energies  of  a  few  hundred  eV  or  less,  which  are  the 
cause  of  Type  A  red  aurora.  In  fact,  the  FI 6  J  sensor  data 
clearly  show  that  besides  a  lowering  of  E  there  is  also  an 
enhancement  in  the  low-energy  electrons  below  0.1  keV. 
This  enhancement  occurs  from  about  1000  UT  until 
2400  UT  encompassing  the  period  of  the  magnetic  cloud. 
The  F 1 5  J  sensor  data  are  also  consistent  with  this  picture  of 
enhanced  low-energy  particle  precipitation  associated  with 
the  magnetic  cloud,  although  these  data  do  not  show  the 
same  enhancement  as  is  seen  in  the  FI 6  orbit. 

[119]  Overall,  the  data  show  the  enhanced  precipitation  of 
low-energy  flux  coincides  with  the  passage  of  the  magnetic 
cloud.  This  is  consistent  with  the  assertion  by  Sivjee  and 
Shen  [1997]  that  interactions  of  magnetic  clouds  with  the 
magnetosphere  may  trigger  the  precipitation  of  low-energy 
particles  which  are  the  cause  of  Type  A  red  aurora.  While 
our  current  observations  provide  considerable  evidence  that 
there  is  a  connection  between  the  passage  of  the  magnetic 
cloud  and  enhanced  low-energy  electron  precipitation,  we 
agree  with  Steele  et  al.  [  1 998]  that  questions  regarding  the 
exact  source  of  such  particles  during  magnetic  cloud  events 
require  a  further  detailed  study  of  the  magnetospheric  data. 
4.2.1.  Shape  of  the  Energy  Distribution 

[120]  Throughout  the  paper,  it  has  been  assumed  that 
Gaussians  are  primarily  responsible  for  the  observed  auroral 
emission.  In  many  previous  studies,  radar  data  had  been 
available  that  showed  that  even  during  superstorms  a 
Gaussian  was  the  predominant  shape  of  the  energy  distri¬ 
bution  [e.g.,  Hecht  et  al .,  1991].  However,  here  no  such  data 
exists  and  so  we  now  show  the  changes  that  occur  in  the 
derived  quantities  if  Maxwellians  were  assumed. 

[121]  Figure  14  shows  plots  of  the  E  and  0/N2  quantities 
assuming  either  a  Maxwellian  with  a  low-energy  tail  or  the 
enhanced  Gaussians  used  above,  both  of  which  are  based  on 
Strickland  et  al.  [1993].  The  left  column  shows  the  time 
series  plots  for  these  derived  quantities  for  two  nights 
19  November  and  20  November.  The  right  column  shows 
the  percent  differences.  On  19  November  (Figure  14a- 14d) 
the  results  show  that  on  average  the  derived  E  and  0/N2 
quantities  are  higher  on  average  by  10-15  and  20-30%, 
respectively.  Since  none  of  the  results  of  this  paper  are 
affected  by  such  an  uncertainty  in  E ,  this  suggests  that 
the  exact  nature  of  the  distribution  is  not  so  significant. 


However,  the  same  cannot  be  said  for  0/N2.  While 
on  17  November  the  essential  conclusions  that  are 
based  on  0/N2  are  not  affected  by  this  uncertainty,  on 

19  November  the  relatively  small  magnitude  of  the  change 
during  the  Joule  heating  event  could  be  affected.  This  is 
especially  so  as  one  could  argue  that  there  is  no  decrease  in 
0/N2,  but  only  a  change  in  the  shape  of  the  distribution. 

[122]  Figures  14e-  14h  show  the  results  for  20  November. 
Here  it  is  seen  that,  although  there  are  periods  where  the 
changes  may  be  larger  (or  smaller)  than  the  nominal 
numbers  quoted  above,  they  would  not  affect  our  conclu¬ 
sions.  If  Maxwellians  occurred,  say,  prior  to  1000  UT,  this 
would  result  in  a  larger  E  and  0/N2.  The  larger  0/N2  would 
thus  indicate  an  even  larger  decrease  than  predicted  by  the 
model  over  the  night.  After  1000  UT  the  choice  of  distri¬ 
bution  has  much  less  of  an  effect  rendering  the  choice  of 
distribution  unimportant  as  to  which  is  used  in  the  analysis. 

[123]  Nevertheless,  there  is  a  method  in  which  the  optical 
data  can  be  used  to  determine  whether  the  auroral  emissions 
are  being  caused  by  Gaussians  or  Maxwellians.  It  had  long 
been  noted  that  Q  and  E  were  often  positively  correlated. 
Christensen  et  al  [1987]  first  showed  that  this  was  a 
consequence  of  the  presence  of  Gaussian  electron  distribu¬ 
tions.  They  showed  that  the  relationship  between  Q  and  E 
should  be  related  by  the  following  when  the  electrons  are 
accelerated  down  a  field  line  producing  a  Gaussians  energy 
distribution 

Q  =  kEg  (1) 

where  g  should  typically  be  between  1.5  and  2  for  most 
acceleration  conditions  depending  mainly  on  the  magnetic 
field  strengths  in  the  magnetosphere  and  ionosphere 
[Christensen  et  al ,  1987].  The  constant  k  is  linearly  related 
to  the  product  of  the  plasma  density  in  the  plasma  source 
region  and  the  square  root  of  the  thermal  energy.  Strickland 
et  al  [1994]  validated  this  using  radar  data  and  showed  that 
Gaussian  distributions  indeed  showed  such  a  correlation 
although  the  g  values  might  be  somewhat  lower  and  higher 
probably  because  the  plasma  conditions  were  changing 
during  the  observation  period.  When  Maxwellians  were 
present  no  correlations  were  seen.  Thus  this  technique 
allows  the  determination  of  not  only  the  type  of  energy 
distribution  but  also  allows  some  insight  into  the  conditions 
in  the  plasma  source  region. 

[124]  Figure  15  and  Table  3  shows  the  results  for  three  of 
the  nights  in  this  study,  17  November,  19  November,  and 

20  November.  The  left  column  of  Figure  15  shows  the  time 
series  plots  of  Q  and  E.  Since  the  E  values  are  not  affected 
much  by  any  prior  knowledge  of  the  distribution,  the 
Gaussian  results  were  used.  It  is  easy  to  see  visually  that 
for  most  periods  the  Q  and  E  values  appear  correlated.  To 
quantify  this,  the  right  column  shows  the  scatter  plots  of  E 
versus  Q,  while  Table  3  shows  the  results  for  k  and  g.  For 
each  night  at  least  one  blue  line  is  shown  that  shows  the  fit 
obtained  by  using  the  fitted  k  and  g  parameters.  Most  of  the 
data  do  appear  correlated  with  g  values  consistent  with  what 
was  found  by  Strickland  et  al  [1994]  confirming  that 
Gaussian  distributions  are  appropriate.  The  following  is  a 
brief  discussion  of  each  night. 

[125]  On  17  November  the  data  from  the  two  substorms 
were  combined  and  shown  as  the  black  data  in  Figure  1 5b. 
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Figure  14.  (left)  Plots  of  Q  and  E  derived  from  red/blue  versus  844.6  nm/blue  technique  at  Poker  Flat 
for  19  November  and  20  November.  These  are  derived  using  the  enhanced  Gaussians  (solid  black 
line)  and  Maxwellians  (dashed  green  line).  Also  shown  are  (right)  percent  differences.  Shown  are 
(a)  19  November  E,  (b)  19  November  percent  difference  for  £,  (c)  19  November  0/N2,  (d)  19  November 
percent  difference  for  0/N2,  (e)  20  November  E ,  (f)  20  November  percent  difference  for  £, 
(g)  20  November  0/N2,  and  (h)  20  November  percent  difference  for  0/N2. 


The  g  value  is  quite  similar  to  the  Strickland  et  al.  [1994] 
results.  They  did  not  report  k  values  though.  The  green  and 
red  data  are  for  the  period  between  the  substorms.  These  do 
appear  appear  to  be  well  correlated.  Thus  Maxwellians 
might  be  appropriate  here.  However,  on  the  basis  of 
Figure  5  this  would  not  change  any  of  the  conclusions  we 
find  with  respect  to  the  observed  versus  modelled  0/N2 
variations. 

[126]  On  19  November  where  the  observed  variations  in 
0/N2  were  small  Gaussians  are  appropriate  during  the 
observing  period  where  the  0/N2  variations  occurred. 
Maxwelllians  might  be  appropriate  at  the  end  of  the 
observing  period.  The  latter  is  where  postmidnight  diffuse 
aurora  can  occur  [Hecht  et  al .,  1995].  Again  a  25%  increase 
in  0/N2  during  the  period  from  1330  to  1500  UT  would  not 
affect  our  conclusions. 

[127]  The  results  from  20  November,  though,  are  the 
most  interesting.  The  analysis  shown  in  Figure  14  indicated 
that  a  knowledge  of  the  type  of  energy  distribution  in  the 
period  before  1 000  UT  would  have  a  much  larger  affect  on 


the  0/N2  results  than  during  later  times.  In  fact  almost  all 
the  data  on  this  night  appear  correlated  and  consistent  with 
Gaussians.  Note  that  during  the  great  magnetic  storm  of 
1986  Gaussians  were  almost  always  observed  [Hecht  et  al ., 
1991]. 

[128]  Unlike  the  other  nights  where  the  scatterplots  show 
the  data  are  clustered,  on  20  November  the  data  undergo 
many  transitions  as  shown  in  Figure  15f  and  in  Table  3. 
This  is  reflected  by  variations  of  g,  which  are  associated 
with  magnetic  field  strengths  in  the  magentosphere  and 
ionosphere,  and,  more  interestingly,  by  variations  in  k 
which  are  related  to  changes  in  the  plasma  source  region. 
From  0600  to  0800  UT  before  the  ICME  shock  hits  the 
magnetosphere,  k  is  quite  low.  These  data  are  shown  in 
black  and  the  blue  fit  is  plotted  through  this  data  set.  After 
the  shock  hits,  another  cluster  of  black  data  are  shown  but 
these  are  shifted  to  the  left.  The  slope  is  quite  similar  but  k 
has  increased  by  somewhat  less  than  a  factor  of  2.  This 
transition  occurs  just  after  the  shock.  The  next  data  set, 
shown  in  green,  represents  the  period  from  930  to  1030  UT, 
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Figure  15.  (left)  Plots  of  the  Blue  intensity  in  kR  (green  line)  and  E  in  keV  (black  line)  derived  from 
red/blue  versus  01(844.6  nm)/blue  photometric  technique  at  Poker  Flat  for  17  November,  19  November, 
and  20  November.  The  Q  and  E  values  are  the  same  as  are  shown  in  Figures  6a,  6c,  8a,  and  8c.  These  are 
derived  using  the  enhanced  Gaussians.  Also  shown  are  (right)  these  E  versus  Q  scatterplots.  In  the 
scatterplots,  different  time  intervals  are  represented  by  different  colors  mostly  indicated  in  the  title.  In 
each  plot  at  least  one  blue  line  is  shown  which  is  the  best  linear  fit  to  the  selected  color  group.  Shown  are 
(a)  17  November  Q  and  E  plots  and  (b)  17  November  scatterplot,  with  black  representing  0800  to 
1 130  and  1330  to  1600  UT,  green  representing  0600  to  0800  UT,  and  red  representing  1 130  to  1330  UT. 
Blue  line  is  a  fit  to  the  black  data.  Also  shown  are  (c)  19  November  Q  and  E  plots  and  (d)  19  November 
scatterplot,  with  black  representing  0900  to  1330  UTand  red  representing  1330  to  1500  UT.  Blue  line  is  a 
fit  to  the  black  data.  Finally,  shown  are  (e)  20  November  Q  and  E  plots  and  (f)  20  November  scatter  plot, 
with  black  representing  0600  to  0800  UT  (Blue  line  is  a  fit  to  the  black  data.),  black  data  to  the  upper  left 
representing  0800  to  0930  UT,  green  data  representing  0930  to  1030  UT,  blue  data  representing  1030  to 
1330  UT,  red  representing  1330  to  1515  UT,  and  light  green  data  representing  1515  to  1600  UT. 


the  transition  period  through  which  the  magnetic  cloud  is 
encountering  the  Earth’s  magnetosphere.  During  this  period, 
the  k  value  undergoes  a  large  increase.  However,  these  data 
still  appear  to  be  connected  to  the  0800-0930  UT  data  and 
the  variations  in  k  and  g  could  be  due  to  changing 
conditions  in  the  plasma  source  region.  The  subsequent 
data  after  1030  UT,  shown  in  blue,  red,  and  light  green, 
clearly  have  a  different  character  than  any  of  the  data  seen 
on  the  previous  nights.  If  these  data  are  simply  interpreted 
in  terms  of  the  Christensen  et  al.  [1987]  analysis  then  the  k 
values  have  increased  by  almost  two  orders  of  magnitude. 


Perhaps  other  processes  are  occurring  since,  as  Christensen 
et  al.  [1987]  note,  their  analysis  does  not  apply  to  down¬ 
ward  field-aligned  currents.  Nevertheless,  these  results 
clearly  support  our  thesis  that  the  magnetic  cloud  period 
is  distinctly  different  than  the  periods  of  auroral  activity  that 
occurred  on  previous  nights. 

5.  Conclusions 

[>29]  In  the  introduction  five  questions  were  introduced 
whose  answers  are  now  presented: 
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Table  3.  Results  of  Fit  to  Function  Q  =  kE*  for  November  2003  Storm  Period 


Day 

UT  Hour  Range 

k 

g 

Plot  Color 

Notes 

17 

0800-  1 130,  1330-1600 

0.14 

1.36 

black 

two  substorms 

19 

0600-0800,  1130-1330 

0.16 

1.05 

black 

two  substorms 

20 

0600-0800 

0.022 

2.14 

black(blue  line  fit) 

prior  to  ICME  shock 

20 

0800-0930 

0.036 

2.23 

black 

after  shock  before  magnetic  cloud 

20 

0930-1030 

0.22 

1.32 

dark  green 

transition  to  magnetic  cloud 

20 

1030-1300 

0.27 

1.77 

blue 

magnetic  cloud 

20 

1300-1515 

0.56 

1.79 

red 

magnetic  cloud 

20 

1515-1600 

1.61 

1.18 

light  green(below  2  keV) 

magnetic  cloud 

[130]  1.  Are  the  energy  distributions  of  the  electron 
precipitation  different  during  the  superstorm  that  began 
about  0800  UT  on  20  November  and  continued  through 
the  first  few  hours  of  21  November,  than  before  the  storm? 

[  13 1  ]  Yes,  although  there  are  two  caveats.  First,  the 
ground-based  data  at  Poker  indicate,  on  the  basis  of  the 
good  correlation  between  Q  and  E ,  that  most  of  the  time  an 
enhanced  Gaussian  distribution  was  appropriate,  rather  than 
a  Maxwellian,  to  represent  the  electron  energy  distribution 
for  the  model  analysis.  Second,  prior  to  1000  UT  on 
20  November  the  E  values  seen  from  both  the  satellite 
and  ground-based  data  were  consistent  with  those  seen  on 
most  prior  days  including  those  seen  on  19  November,  the 
geomagnetically  quietest  day  of  this  period.  However,  after 
1000  UT  on  20  November  and  continuing  until  the  first  few 
hours  of  2 1  November  the  E  values  are  lower  than  during 
other  periods,  including  the  period  of  enhanced  auroral 
activity  seen  on  17  November.  The  low  E  values  are 
obtained  at  the  ground-based  sites  in  Alaska  and  Canada 
and  also  in  most  of  the  DMSP  and  POES  satellite  data. 
While  occasionally  some  of  the  satellite  data  showed  low  E 
values  on  17  November  and  19  November,  the  period  from 
about  1000  to  2400  UT  on  20  November,  unlike  any  other 
period,  had  electron  energy  distributions  characterized  by  a 
large  enhancement  in  the  flux  below  0.1  keV.  The  early 
hours  of  21  November  showed  some  enhancement  in  this 
low-energy  flux  but  much  reduced  from  that  seen  on 
20  November. 

[132]  2.  Is  there  a  connection  between  the  passage  of  the 
magnetic  cloud  during  this  storm  and  the  low-energy 
particle  precipitation  as  suggested  by  Sivjee  and  Shen 
[1997]  from  earlier  observations? 

[133]  Although  the  20  November  superstorm  was  initiated 
by  the  passage  of  a  solar  wind  shock  by  the  Earth  near 
0800  UT,  a  decrease  in  the  average  electron  energy  or  a  large 
enhancement  in  low-energy  particles  was  not  associated  with 
the  shock.  A  large  enhancement  in  the  very  low  energy 
electrons  was  observed  by  the  DMSP  F16  J  sensor  during 
the  period  of  1000  to  2400  UT  on  20  November,  which 
corresponds  to  the  time  of  the  passage  of  the  magnetic  cloud 
as  defined  by  Gopalswamy  et  al.  [2005b].  Some  lesser 
enhancements  compared  to  the  previous  nights  were  seen 
just  before  and  after  the  period  of  the  cloud  suggesting  the 
cloud  boundary  may  not  be  sharp.  In  addition,  during  the 
magnetic  cloud  passage,  large  red/blue  ratios  were  seen 
suggestive  of  the  presence  of  Type  A  red  aurora.  Thus  taken 
together  these  support  the  assertion  by  Sivjee  and  Shen 
[1997]  that  the  interaction  of  the  magnetic  cloud  with  the 
magnetosphere  may  trigger  an  enhancement  of  low-energy 
electron  precipitation  responsible  for  such  aurorae. 


[134]  Nevertheless,  the  exact  connection  between  such 
precipitation,  the  passage  of  the  magnetic  cloud,  and  any 
delayed  effects  of  the  large  shock  that  occurred  about  0800 
UT  need  to  be  further  investigated.  To  date,  the  connection 
between  low-energy  particles  and  magnetic  clouds  rests  on 
only  the  example  in  this  paper  and  the  previous  few  studies 
by  Sivjee  and  Shen  [1997]  and  Steele  et  al  [1998].  In  this 
work  we  have  shown  using  the  measured  relationships 
between  Q  and  E  that  the  period  after  1000  UT  on 
20  November  is  distinctly  different  than  what  was  seen 
prior  to  that  period.  In  the  context  of  the  Christensen  et  al. 
[1987]  analysis  this  might  mean  changes  in  the  plasma 
source  region,  or  perhaps,  different  processes  are  occurring 
than  the  electric  field  acceleration  considered  in  that  work. 

[135]  3.  To  what  extent  do  the  0/N2  model  (NRLMSIS) 
predictions  and  actual  observations  agree  during  the  period 
before  and  during  the  superstorm? 

[  1 36]  Although  NRLMSIS  predicted  changes  in  0/N2 
generally  follow  Kp  because  of  its  3  h  timescale,  the  model 
does  not  fully  capture  composition  changes  due  to  geomag¬ 
netic  activity.  During  the  low  geomagnetic  activity  night 
of  19  November  the  model  and  data  results  are  similar 
reflecting  the  small  changes  that  occurred  that  night.  On  the 
superstorm  night  of  20  November  the  model  captures  the 
general  trend  of  decreasing  0/N2  during  the  night  but 
underpredicts  the  magnitude.  The  data  show  a  factor  of 
three  decrease  in  0/N2,  which  is  comparable  to  what  was 
found  during  the  1 986  superstorm  discussed  by  Hecht  et  al 
[1991].  However,  the  largest  difference  is  on  17  November 
where  the  model  and  data  show  large  differences  during  the 
time  of  peak  fluence  Fy.  This  difference  is  most  likely 
related  to  the  shorter  time  duration  of  the  Joule  heating  on 
this  night  compared  to  the  sustained  periods  of  intense  Joule 
heating  on  20  November. 

[137]  4.  In  particular,  do  we  see  any  difference  in  0/N2 
between  the  model  and  observations  during  the  time  period 
of  the  magnetic  cloud  event? 

[1 38]  No.  The  magnetic  cloud  period  does  not  show  any 
particular  increase  in  the  rate  of  decline  of  0/N2.  The  only 
interesting  aspect  is  that  the  dayside  decrease  in  0/N2 
reported  by  GUVI  appears  to  begin  at  the  time  of  the 
passage  of  the  magnetic  cloud. 

[139]  5.  Do  any  of  the  observed  atmospheric  composition 
changes  correspond  to  inputs  of  auroral  Joule  or  particle 
energy  into  the  thermosphere? 

[140]  Yes.  The  decreases  in  0/N2  appear  to  scale  with  the 
energy  input  into  the  atmosphere  over  periods  of  an  hour  or 
more.  The  largest  decreases  in  0/N2  occur  during  periods 
when  Fy  peaks.  Fy  especially  exceeds  the  fluence  due  to 
auroral  particles  on  20  November.  This  is  similar  to  what 
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has  been  previously  reported  and  is  summarized  by  Hecht  et 
al  [2006]. 
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